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ABSTRACT 

We present results for a galaxy-galaxy lensing study based on imaging data from 
the Canada-France-Hawaii Telescope Legacy Survey Wide. From a 12 million object 
multi-colour catalogue for 124 deg^ of photometric data in the u*g'r'i' z' filters, we 
compute photometric rcdshifts (with a scatter of crAz/(i+z) = 0.033 and an outlier rate 
of 77 = 2.0 per cent for i' ^ 22.5) and extract galaxy shapes down to i' = 24.0. We 
select a sample of lenses and sources with 0.05 < Zd ^ 1 and 0.05 < Zg ^ 2. Wc fit 
three different galaxy halo profiles to the Icnsing signal, a singular isothermal sphere 
(SIS), a truncated isothermal sphere (BBS) and a universal density profile (NFW). 
We derive velocity dispersions by fitting an SIS out to 100 kpc to the excess sur- 
face mass density AE and perform maximum likelihood analyses out to a maximum 
scale of 2 Mpc to obtain halo parameters and scaling relations. We find lumi- 
nosity scaling relations of di-od iO-24±o.03 ^Y\c red lens sample, iTbiuo c>c iO-23±o.03 
for blue lenses and a oc 7^o.29±o.02 ^-^^q combined lens sample with zeropoints of 
a*^^ = 162 ± 2 km s~\ a^^^^ = 115 ± 3 km s'^ and a* = 135 ± 2 km s'^ at a chosen 
reference luminosity L*, = 1.6 x lO^'' h^^ Lr\Q- The steeper slope for the combined 
sample is due to the different zeropoints of the blue and red lenses and the fact 
that blue lenses dominate at low luminosities and red lenses at high luminosities. 
The mean effective redshifts for the lens samples arc (zrod) = 0.28 for red lenses, 
{zh\uc) = 0.35 for blue lenses and (z) = 0.34 for the combined lens sample. The BBS 
maximum likelihood analysis yields for the combined sample a velocity dispersion of 
(7* = 131^2 km s~^ and a truncation radius of s* = 184^^4 h^^ kpc, corresponding 
to a total mass of M^^^^j ggg = 2.32to:25 x 10^^ ^"^ ^^o and a mass-to-light (M/L) 
ratio of Mj*^^^^) ggg/L* = 178l^g /i Mq/L^',© at L*,. At a given luminosity, both 
velocity dispersion a and truncation radius s are larger for red galaxies than 
for blue galaxies. For an NFW profile, we measure at L*, a virial radius of 
''loo ~ I33I2 kp*^ ^ concentration parameter of c* = 6.4^q y, implying a 
virial mass of Afjoo = '''•S^o's x 10^^ h^^ Mq. At L*, for blue galaxies the concentra- 
tion parameter is slightly higher than for red galaxies and r2oo is significantly lower. 
For the combined sample, if described as a single power law, the M/L-ratio scales 

as A/total, BBs/i OC L the concentration parameter scales as c cx L -o n. 

Analysing the M/L-scaling for red and blue galaxies separately, we find that a broken 
power law (with a flat slope at high luminosities) provides a more appropriate de- 
scription for the red and possibly also for the blue galaxies. We measure M2oo/Afstar 
for red galaxies over 2.5 decades in stellar mass. We find a minimum for this ratio at 
Mitar ~ 3 — 4 X 10^° Mq with a strong increase for lower stellar masses. 

Key words: gravitational Icnsing: weak - galaxies: distances and redshifts - galaxies: 
haloes - galaxies: kinematics and dynamics - dark matter 
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1 INTRODUCTION 

Dark matter is the dominant mass component in the Uni- 
verse, not only on scales of clusters or groups of galaxies but 
also on galaxy scales where dark matter haloes surround the 
luminous baryonic cores of galaxies. This brings up the ques- 
tion of how to measure and quantify the properties of dark 
matter galaxy haloes. On short scales the luminous matter 
can serve as a tracer to investigate the mass distribution 
with dynamical methods (e.g. measuring rotation curves of 
spiral galaxies or velocity dispersions of early type galax- 
ies). On lar ger scales one can us e satellite galaxies as trac- 
ers (see e.g. IPrada et al.ir2003l or lConrov et al.ir2005l ) or can 
study the gravitational lens effect that dark matter haloes 
impose on background galaxies: they introduce a small co- 
herent distortion to the shapes of background galaxies with 
an amplitude depending on the angular diameter distances 
and relative positions of foreground and background galax- 
ies and on the dark matter distribution of the galaxy haloes. 
As the gravitational impact of a single galaxy is too 
small to be detected individually, a large sample of 
foreground and background galaxies is required to mea- 
sure this effect statistically. The first successful detec- 
tion of galaxy-galaxy-lensin g (GG L) was reported by 
iBrainerd, Blandford &: Smaill (|l996l ). although the analy- 
sis was limited by t he small area covered by th e obser- 
vations, followed by Ideir Antonio fc TvsonI (|l996l ) in the 
Hubble Deep Field (HDF). Later on GGL measurements 
from space wit h the Hubble S pace Telescope (HST) have 



been done b y [Griffiths et all (11 996. HST Medium Deep 



Survey) and iHudson et al.l (|l99a . HDF). Fu rther ground- 
based GGL analyses have been c arried out by Wilson et al.l 
l|200lll with CFHT/UH8K or bv lHoekstra et al.l (|2003h us- 
ing data from the Canadian Network for Observational 
Co smology Field Galaxy Re d shift Survey (CN0C2) and 
bv iHoekstra. Yee fc Gladderd (|2004l ) investigating a signif- 
icantly larger area of 45.5 deg^ Red-Sequence Cluster Sur- 
vey (RCS) imaging data. While in bot h cases photometric 
redshift information was not available, iKleinheinrich et al.l 
lj2006h analysed Wide Field Imager (WFI) data from 
the 'Classifying Objects by Medium-Band Observations 
in 17 filters' (COMBO-17) survey, measuring th e scal - 
ing for the TuUy-Fisher relation (jTuUv fc Fisheil Il977l ) 
and investigating the halo properties assuming an NFW 
profile. Up to today the largest survey area used to 
quantify dark matter properties on gala xy scale i s the 

200d . 

or IScliak '200d). ISheldon et al] 



Sloan Digital Sky Sur vey (S DSS, e., 
McKav et all I2001I ' 



Mandelbaum et al 



2004h 



2006g): Mandelbaum. Seliak fc Hiratal 



l|2008h and ISchulz. Mandelbaum fc PadmanabhanI (120101 ) 

measured the excess surface density AE and mass profiles 
for various samples using spectroscopic redshift informa- 
tion for the lenses and a combination of photometric red- 
shift and stat istical redshift distribut ions for background ob- 
jects. Finallv. fvan Uitert et al.l (j201lh combined high quality 
imaging data for the background shape estimation from the 
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RCS2-r-'-band data and the profound knowledge of the fore- 
ground redshifts available from the SDSS DR7 to quantify 
the dark matter halo properties as a function of luminosities, 
stellar and dynamical masses. 

In this study we make use of 124 deg^ pubhc CFHTLS- 
Wide imaging data, covering u*g'r'i'z' photometry, thus 
providing an excellent base for photometric redshift estima- 
tion and perform a GGL analysis. We give an introduction 
to the theory of GGL and introduce our notation in Sec- 
tion [2] In Sections 13.11 and 13.21 we discuss the imaging and 
the spectroscopic data used in this paper. The creation of 
colour catalogues and the computation of photometric red- 
shifts are presented in Section 13.31 The background shape 
estimation is described in Section [3.41 We present the prop- 
erties of our sample catalogues in Section 13.51 define lens 
and background source samples in Section 13.61 and discuss 
the lens mass errors introduced by possible errors in shape 
measurement and photometric redshift estimation in Sec- 
tion |XZ1 The weak lensing analysis, including measurement 
of the tangential shear 7t, measurement of the excess sur- 
face mass density AE and maximum likelihood analyses for 
various subsamples concerning spectral energy distributions 
(SEDs), luminosities and environments, is described in Sec- 
tion|3| In SectionOwe summarise and conclude our analysis. 
In the Appendix we measure the systematic errors to verify 
the integrity of our analysis and compare our lensing results 
to simulations of the 3D-line-of-sight(LOS)-projected weak 
lensing signal predicted by our results, i.e. using the mea- 
sured best-fit parameters and luminosity scaling relations 
for the dark matter haloes. 

Throughout this paper we adopt a cosmology with 
f^M = 0.27, = 0.73 and a dimensionless Hubble parame- 
ter h = 0.72, unless explicitly stated otherwise. All referred 
apparent magnitudes and colours are given in AB, all rest 
frame magnitudes are calculated in vega system, assuming 
a Hubble constant of Ho — 72 km s~^ Mpc~^. 

Our created catalogues (photometry, photometric red- 
shifts and shears) are publicly availablelj 



2 THEORY OF GALAXY-GALAXY LENSING 

Weak gravitational lensing describes the impact of line- 
of-sight structures on the shape, size and magnitude of 
background galaxies in the limit where no multiple images 
are produced. If the galaxy shape distortions are small, 
the shape distortions can be directly related to the density 
field producing them. The inversion of the shear-mass 
density relation into the mass density-shear relation is the 
Kaiser and Squires equation ([Kaiser fc Sauirei Il993l ). It 
has paved the way for weak lensing analyses in general (for 
the equations we refer the reader to Section [2. l|l and it also 
states the principle of weak lensing analyses using galaxy 
shapes: The mass density field can be inferred from galaxy 
shapes because they carry the imprint of the tidal field of 
the gravitational potential. 

Weak gravitational lensing analysis techniques can be 
divided into three branches: one can obtain mass maps for 
individual galaxy clusters, one can study the mean radial 
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profile of an ensemble of clusters, groups and galaxies, 
and finally one can measure the statistics of dark matter 
fluctuations as a function of scale. For the first case one 
needs a high signal-to-noise lensing signal, i.e. a massive 
cluster and deep observations, providing a high number 
density of background galajcies. One can then turn a 2D 
gravitational shear map into a 2D mass map without 
requiring a reference to the light distribution of the cluster. 
For the second analysis method the signal of the individual 
structures can be lower. The stacking or averaging tech- 
nique then yields the mean properties of the sample. The 
stacking requires a reference position, which is taken to 
be, e.g. the brightest cluster galaxy in case of cluster and 
group weak lensing and the centre of light for the galaxy 
of interest in case of GGL. The analysis therefore then 
results in a correlation between galaxy positions and the 
mass associated to them (galaxy mass correlation function). 
The third method (cosmic shear) analyses the correlations 
of gravitational shears along pairs of light rays traced by 
galaxies. The shear-shear correlation function has a very 
low amplitude, since it is the product of two small numbers. 
Therefore cosmic shear analyses require to analyse the 
largest amount of galaxy pairs to obtain a significant result. 
If one ignores the problem of measuring galaxy shapes, 
then the two remaining sources of bias in weak lensing 
measurements are due to incorrectly locating the galaxies 
along the line-of-sight and due to limited validity of the 
assumption that galaxies are randomly oriented on the sky. 

In this Section we shortly summarise the theory of 
GGL in Section 12.11 and introduce the tangential alignment 
7t and the excess surface density AE. In Section 12.21 we 
describe the galaxy halo profiles investigated later on. 
Section 12.31 introduces the scaling relations which we have 
assumed (and later on show to hold) when we analyse the 
profile halo parameters. In Section 12.41 we investigate the 
influence of a group or cluster halo on the lensing signal of 
a galaxy. Finally in Section [2.51 we introduce the maximum 
likelihood analysis, which is a more sophisticated method to 
analyse the GGL signal, since it fully accounts for multiple 
deflections of light rays by multiple foreground galaxies. 



2.1 Tangential alignment and lens mass excess 
surface density 

We briefly re view some relatio ns needed later on. We refer 
the reader to ISchneideJ (|2006 ^) for an extended overview. 
In gravitational lensing the light deflection angle a can be 
obtained from a deflection potential ip as a = Vip, and the 
lens equation reads j3 = d — Ct = d — Vi/'- The Jacobian de- 
scribes the mapping of small lig ht bundles (|Sch ncidcr 200S; 
ISeitz. Schneider fc Ehler3ll994l l and depends on the conver- 
gence 



2 ^ Ec 



and the gravitational shear 



^ 7i + *72 = 2 (^.11 



- ip,22) + iip,: 



(1) 



(2) 



K is the ratio of the surface mass density E of the lens and 
the critical density 



47rG Di L>ds ' 

depending on the angular diameter distances and Ds 
from the observer to the lens and source and Dds from the 
lens to the source. The gravitational shear 7 is a two compo- 
nent quantity describing the size and direction of the tidal 
forces distorting the shapes of light bundles. Equations 
and ((2]) show that k and 7 are rela t ed to each other via the 
deflection potential. In fact iKaiseil l|l995h has shown that 



71,1 + 72,2 
72,1 — 71,2 



(4) 



holds and that one can obtain Vk from shape esti- 
mates in the very weak lensing limit, k <C 1. For 
the more general case relati o ns are given in iKaiseij 
(j 19951 ), [Sciineider fc Seitd (Il995l). [Seitz fc Schneided (ll995D 
and are used e.g. in iKaiser. Squires fc BroadhursTlll995h , 
ISeitz fc SchneideJ l|l996l ) and ISeitz fc Schneidej l|200li rto 
obtain K-maps from local shape estimates. 
We will now focus on the information that can be extracted 
from GGL, where the signal of many galaxies is stacked and 
only the mean t ange ntial shear pro flle can be observed. 
ISchneideil l|200^ and'K aisej (|l995l ) have shown that for an 
arbitrary mass distribution the mean tangential shear (jt) 
on a circle at radius 8 can be written as a function of the 
mean surface mass density R within that circle and the mean 
surface mass density (k) on the edge of that circle: 



(70 (0) = k(0) - (n) (6) 



(5) 



For circular mass distri butions this r elation was derived by 
iBartelmannI (|l995l ) and lKaised (|l996l ). 
Equation ([5]) can be rewritten as 



Ec {7t) (R) = S(i?) - (E) (R) = AE(i?) 



(6) 



which is called the excess surface mass density AE. This 
means that if for a lens system the source and lens dis- 
tances and the m ean tangential shear is known then AE(ii) 
can be derived. iMcKav et all l|200ll ) were the flrst who 
made use of this relation in the context of GGL. They di- 
rectly measured AE for galaxies as a function of the dis- 
tance R to the galaxy centres instead of just presenting the 
tangential gravitational shear (7t) as a function of angu- 
lar scale 0. Before that the mass density associated with 
galaxies was derived from the observations using paramet- 
ric models for the gravitational shear and statistical de- 
scriptio ns for the dis tribution of background and fore ground 
objects terainerd. Blaiidford fc Srnail 1996; Hoekstraet al.l 
l2003l : iHoekstra. Yee fc Gladder j|2004l : iParker et aLlbOol v 

2.2 Profiles for galaxy haloes 

The measurements of AE(7?) (equation |6]) relative to the 
centres of galaxy light directly link galaxy evolution to struc- 
ture formation. In the most simple view a galaxy can be 
pictured as being composed of an ensemble of stars in the 
centre surrounded by a dark matter halo which is described 
by an analytic profile. Since the dark matter halo is sup- 
posed to belong to the galaxy it has to be finite in extent 
and in mass. 
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If galaxy haloes were all isolated one could just measure their 
profile directly, in particular how haloes end or are cut off in 
the outskirts. In a more complete picture one however has to 
take into account that galaxy haloes might be embedded as 
satellites in more massive galaxy, group or cluster haloes, or 
that their halo could contain satellite galaxies itself. There- 
fore, especially at larger distances, the individual profiles 
barely can be traced as we additionally measure the parent 
haloes in which the investigated galaxies are embedded in, or 
we measure the line-of-sight projection effects of other haloes 
within the same structure (e.g. in group and cluster cores). 
Furthermore, multiple gravitational deflections by different 
line-of-sight structures introduce further difficulties, which 
need to be tak en into account in a quantitive analysis (see 
lBrainerdll2O10l '). 

Although the situation seems complex it makes sense to 
nevertheless use simple halo profiles and constrain their pa- 
rameters: On scales smaller than 200 kpc, the galaxy 
mass correlation function is dominated by the halo in whi ch 
the galaxy resides (see Sec t ion [2:41 and ISeliak et al.ll2005l or 
iMandelbaum et al.|[2005bl . fig. 1). This can be used to es- 
tablish scaling relations and thus to model simultaneously 
galaxy haloes with various luminosity reasonably well. Then 
one can quantitatively describe all haloes that host a galaxy 
in their centre and vary their joint parameters such to match 
observations best. In this way only haloes which do not host 
a galaxy or which host a galaxy fainter than the detection 
limit get ignored. In the following we will model each galaxy 
halo as a function of the galaxy luminosity and scaling re- 
lations relative to a halo with a fiducial luminosity L* . The 
motivation to analyse the galaxies as a function of luminos- 
ity and to assume and to investigate scaling relations are the 
following: 

(i) We do not only want to analyse 'the' average dark 
matter halo (a result we could easily obtain by just stacking 
the signal with respect to all foreground galaxies). 

(ii) We want to establish the link between dark matter 
halo properties and the optical properties of galaxies. The 
most basic galaxy property is the luminosity. In the redder 
bands it can also be taken as a proxy to stellar mass, which 
then allows to study how effectively baryons are turned into 
stars as a function of dark matter halo depth. 

(iii) The galaxy luminosity is tightly related to the total 
halo potential at least for the central regions of a galaxy, 
partly because the stars are a major contributor to the po- 
tential in the centres of a galaxy, partly because of the lumi- 
nous and dark matter conspiracy. This is causing the transi- 
tion in the haloes from luminous to dark matter dominance 
to be such that spiral rotation curves are fiat and ellipticals 
have a close to isothermal mass profile. From this picture we 
(in agr eement with the Tully- Fish er and Faber- Jackson re- 
lation, |^berTjackso3[l97i and lTullv fc Fisheij[l977l ) ex- 
pect the halo parameters (in particular the halo 'depth') to 
scale with luminosity. 

(iv) If the halo parameters scale with luminosity we can 
assume scaling laws when analysing the GGL signal and 
obtain the halo parameters for a reference halo with high 
signal-to-noise by jointly analysing galaxies with various lu- 
minosities. 

(v) We can also investigate (by splitting into luminosity 
subsamples) whether the assumed scaling laws indeed hold 



and measure e.g. the coefficients for the power law scaling 
relations between luminosity and halo parameter. 

We will not distinguish between central galaxies and satel- 
lites at this point. This means that our results will de- 
scribe the mean dark matter halo of a galaxy as a func- 
tion of its luminosity. Then a galaxy which is in the cen- 
tre of a halo is expected to have a larger mass than 
our mean, and a galaxy which is a satellite is expected 
to have a lower mass (due to halo stripping) than this 
mean. The halo profiles investigated are the singular isother- 
mal sphere (SIS), the truncated is othermals sphere (BBS, 
iBrainerd. Blandford fc Smaillll996l) and the universal dark 
matter profile of NFW l|Navarro, Frenk fc Whitelll996l ). 



2.2.1 Singular isothermal sphere (SIS) 
The singular isothermal sphere (SIS) 

Psis(-R) 



(7) 



(8) 



27rGi?2 

has a projected surface density of 

^2 

^^'^^^^ = 2GR ' 

with a being the line-of-sight velocity dispersion of test par- 
ticles (stars) and R the distance from the centre of mass 
in the plane of projection. For an SIS the convergence and 
tangential shear are equal. 



7t,sis(6')=Ksis(6')=47r(^) 



261 



(9) 



where we have used R — 9 D^. The line-of-sight projected 
mass within a sphere of radius R is given by 



Msis(< R) = ^R , 
which means that the total SIS mass diverges. 



(10) 



2.2.2 Truncated isothermal sphere (BBS) 

'Brainerd. Blandford &: Smaill (| 19961 ) introduced the trun- 
cated isothermal sphere (BBS) with a volume density of 



PBBS{R) 



(11) 



which has one further free parameter compared to the SIS, 
the truncation radius s. For infinite truncation radii s the 
density profile converges to that of an SIS. On short scales 
[R <C s) the difference to an SIS is negligible but at a radius 
of R = s the volume density has dropped down to only half 
the value of an SIS with the same velocity dispersion. For 
larger distances R^ s the BBS volume density profile falls 
off oc _R~*. The surface mass density equals 



SBBs(i?) = 



2GR 



R 



VR^ 



(12) 



which can be rewritten with 9^ — 9/s in terms of the dimen- 
sionless convergence kbbs as 



HBBS{9) = 
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V9^TW 



(13) 
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The gravitational shear equals 

^E.sis f + 29s 6'^ + 26'i 



7BBS {d) 



29 



(14) 



The line-of-sight projected mass within a sphere of radius R 
is given by 



A2bbs(< R) = „ arctan — 
G V s 



adding up to a total halo mass of 

TTa^s 

Altotal.BBS 

7.3 X 10'^ /i"' Mq 



(15) 



G 
a 



lOOOkms-i/ VlMpc 



(16) 



2.2.3 Universal density profile (NFW) 

The third profile considered is the universal density 
profile, also known a s NFW pr ofile, introduced by 
iNavarro. Frenk fc White! j 19961 , Il997l ) with a volume den- 
sity of 



Pnfw(^) = 



ScPc 



{R/r,){l + R/r,Y 



where 



Pc 



?,H{zf 

8-kG 



(17) 



(18) 



stands for the critical density of the Universe at redshift z. 
The NFW profile has two parameters, the scale radius rs, 
which can be expressed in terms of the radius where the den- 
sity profile turns from /9nfw(^?) oc R~^ to pNFw(^?) oc R~^ 
and the density contrast 



200 



3 ln(l -f c) - c/(l + c) ■ 
c is called the concentration parameter and defined as 

^200 



(19) 



(20) 



giving the ratio between the virial radius r2oo (within which 
the mass density equals 200pc) and the scale radius rs. 

The projected surfa ce mass dens i ty of an NFW pro- 
file has been derived in iBartelinann lll996|), the g r avita - 
tional shear can be inferred from Wright fc Brainerdl (|200G| ) 
or ca n be calculated str aightforwardly with equations 7-11 
from iBartelmannI (|l996l ). Adopting a dimensionless radial 
distance x = R/ra = 9/9i^ we can write the convergence as 



Ec(a:^-l7 
2r3fcPc 
3 

2rs<5c 



=artanh 



_ Eo(a;2_l) 



= arctan 




The gravitational shear is given by 



+4 1n(i)] 
--i^g>{x) 



X <1 
x = 1 

X <1 



(21) 



(22) 



The functions g<{x) and g>{x) are independent from cos- 
mology and do only depend on the dimensionless radial dis- 
tance x: 



9<(x) = 



artanh^ (1 



(a;2 - 1) 



+ 



4 artanhY/(l — x)/{l + x) 
(a:2 - 1)(1 -3;2)i/2 



(23) 



arct an ^/{x — + x) 4 



x^\/x^ — 1 

2 ^ 4 arctanV(T - 1)/(1 x) ^^^^ 



(x^-l) ' (a:2_ 1)3/2 

The virial mass, the mass within the r2oo, is given by 

_ , , SOOtt 3 

M200 = AiNFw(< ^200) = ~~^P'' "^200- 

The rotation velocity at the virial radius is 



^200 



GM200 
'"200 



SOOvrpcG 



'"200 



(25) 



(26) 



If we fit BBS and NFW profiles to a mass distribution, 
we do not expect that the total BBS- and the Af2oo-NFW- 
masses are equal: the M2oo-NFW-mass only gives that part 
of the halo mass which is already virialised, and does not 
specify the mass associated with the structure outside this 
radius. It however makes no sense to integrate this profile 
for much larger radii, since th e total mass would diverge. 
iBaltz. Marshall fc Oeurl (|2009l ) presented a modification of 
the dark matter profile that is NFW-like inside the virial 
radius and describes a cutoff relative to the original NFW 
profile for larger radii. The gravitational shear then falls off 
faster than for the NFW-profile (see their fig. 2). To mea- 
sure these truncation details is not the goal of this work. We 
however point out that these truncated NFW haloes are also 
described approximately by a BBS halo. The total mass of 
the truncated NFW halo then is similar in interpretation to 
the total BBS mass, i.e. it equals the total mass associated 
with the halo. Therefore we expect that M2oo-NFW-masses 
are lower than BBS masses obtained from our model fits. 



2.3 Scaling relations for the galaxy halo profiles 

We assume the following scaling relations between halo pa- 
rameters and galaxy luminosity: 



L \ 1M 



/ M/L 
\{M/LY 



iw) - (f) 

L \Vm/l 

JT ) ' Vm/l 



riM 



(27) 



(28) 



In these two upper equations M can be either the virial mass 
M200 of the NFW or the total mass Mtotai,BBs of the BBS 
profile. The corresponding scaling indices »7a/2oo ^-^d Wbbs 
can be different, in principle. The luminosity in this relation 
is not the bolometric one but that for a specific filter, which 
we choose to be the r'-band filter in the following. 
The halo velocity dispersion vs. luminosi ty relation is pa- 
rameterised by the Faber-Jackson relation (|Faber fc Jackson! 
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I TuUv fc Fisiie3ll97'i1 ) 



1 19761) or (using tiro t oc a) the TuUy-Fisher relation 

(29) 



(30) 



.cr'/ \L* 

The halo-size vs. luminosity relation is described by 

L \vs 



is") (l* 



Inserting equations ([27}, p9p and (|30p into equation ([161 
we obtain 



??A/bbS = "^Vt + Vs 



(31) 



Using the SDSS-g'-band for the mass-luminosity rela- 
tion and analysing the GGL- signal in an early SDSS- 
data set iGuzik fc SeliakI ([20021 ') have found ■r]M = 1.2 ± 
0.2, or riM/L = 0.2. This is also the scaUng behavior 
found for the dynamical mass-to-light ratio for the cen- 
tres of ell iptical galaxies, the so-cal l ed fu nda mental plane 
(see e.g. [Bender. Burstein &: Fabeil Il992l or [Saglia et al.l 
I2OIOI '). The TuUy-Fisher and Faber- Jackson power law 
indices that have been measured in various filters and 



galaxy samples (see e.g. [Payies et al. 19831 ana. l ysing 
faint early type galaxies and iMatkovic fc GuzmanI |2007[ 
an alysing dwarf early type galaxies, both in _B-band, 
or [Nigoche-Netro et al.l I2OI0I analysing SDSS early type 
galaxies in g- and r-band for the Faber- Jackson relation 
and analysing the Tully-Fisher relation for disc galaxies 



e.g. iBamford. Aragon-Salamanca Sz Milvang- Jensen! 2006 



Miller et al.l I2OIII in _B-band and [Pizagno et al.l 12005 



Fernandez Lorenzo et al] [20091 or iReves et al.l 201ll in 
7-band). The 'classica l' Faber- Jackson index is 0.25 
(|Faber fc Jackson! 1 19761 ). but variations depending on the 
sample selection, e.g. luminosity interval and galaxy mix- 
ture ('pure ellipticals', SOs, spirals or a combination) have 
been reporte d (see e .g. Williams. Bureau fc Cappellarill2010l 
or [Cappclla ri et al.l |2012| ). When Equation (|29p was as- 
sumed or attempted to be measured in the analysis of 
the GGL signal different assu mption s were made and dif- 
ferent results were obtained. ISeliakI ([2002D . for instance, 
followed the classical Faber - Jackson relation for elliptical 
galaxies, iHoekstra. Yee fc Gl adders (2004) assumed a scal- 
ing of fj oc L° '^ for a mixed SED sample. On the other hand 
iKleinheinrich et al.l (|2006l ) measured values of rja- between 
0.3 and 0.4, depending on the maximum considered pro- 
jected separations between source and lens. 
If one assumes tim = 1.2 and rjcr ~ 0.3 to hold then rjs — 0.6 
must hold as well. These three values provide kind of a 
generic model to analyse a combined galaxy sample and ob- 
tain the velocity dispersion and size of the 'reference' halo 
with luminosity L* . Our analysis will go beyond these as- 
sumptions and aims to measure these parameters for the 
scaling relations. In addition we'll measure the scaling rela- 
tions and the values for the reference halo with luminosity 
L* for the red and blue galaxies separately. 
We write the virial radius vs. luminosity relation of the NFW 
profile as 



/ r2oa \ _ / L \ -n^2 
V r?nn / \L* ) 



This implies 



'7M2 



3?7r: 



(32) 



(33) 



The zeropoint of this relations changes with redshift as the 
critical density is redshift dependent as well (see equa- 
tions [TH| and [^ . 

For the con c entrat ion-mass-relation we take the relation of 
iDuffv et al.l (|2008l '). derived for simulated haloes made of 
dark matter only and which depends on the redshift of the 
halo, to hold, 

,^-0.084±0.006 , N-0.47±0.04 ,nA\ 
C CX M200 iX + z) ■ (34) 

Already e.g. [Bullock et all (|200ll ) or IShaw et all (|2006D 

found a concentration-mass-relation, implying a slight de- 
crease of the concentration parameter with increasing mass. 
In the past years also furt her mass-concentration rel ations 
have been discussed by e.g. iBhattacharva et al.l (|2013l ). who 
found a relation with slightly higher a mplitude, but very - 
si milar scaling behavior as the one o f iDuffv et al.1 (|2008D 
or iKlvpin. Truiillo- Gomez fc Primac^ (|201lf ). reporting an 
even higher amplitude a nd a shallow e r decr ease with in- 
creasing mass. Recently [Prada et al.l (|2012l ) presented a 
concentration-mass-relation, including a novel feature, as 
above a virial mass of approximately lO'^^ Mq the con- 
centration parameter begins to increase again instead of fur- 
ther decrease. 

Using the M2oo-niass and luminosity relation in equa- 
tion (|27p we obtain the concentration-luminosity relation 



with 



nc = 



-0.084 ± 0.006 



(35) 



(36) 



2.4 Influence of a cluster or group halo on AE 

We investigate the influence of a parent cluster, group or 
massive galaxy halo on the observable excess surface mass 
density AE of a galaxy. In this case the galaxy halo is 
a satellite in a more massive halo. This effect has b een 
investigated and il l ustrate d by iGuzik fc SeliakI ([20021 ') or 
iMandelbaum et al.l ([2006d ) in the past already. Since we 
need some numbers for later on, e.g. on which scales this 
effect is important, we give a relevant example here as well. 
We assume that a galaxy halo described by an SIS with 
velocity dispersion a — 120 km s~^ is embedded in a 
group halo described by an NFW profile with concentration 
c = 5 and r2oo ~ 600 kpc. We place this galaxy halo 
centrally in the group halo and with projected distances of 
0, 7, 100, 200 and 400 kpc. The total excess surface 
densities AE with respect to the galaxy halo centre are 
shown as a function of projected distances in Fig. [T] with 
black curves (upper two, middle two and lower left panel). 
The contribution of the galaxy halo and the group halo to 
the combined signals are shown in red (dashed) and green. 
If the galaxy halo is slightly offcentred from the group 
halo, a turnover of the combined excess surface density at 
small radii can be observed. The farther the galaxy from 
the group centre, the larger the scale at which the excess 
surface density of the combined signal equals the excess 
density of the galaxy halo. If one considers an ensemble of 
galaxies that is homogeneously spread over the group halo 
out to radii of 600 kpc (this was assumed to be the 
case in the lower right panel of Fig. [T]) then the group halo 
contribution averages out (i.e. if limited to small enough 
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Figure 1. Observable excess surface mass density AS for galaxies embedded in a group or cluster halo. For the galaxy halo we use an 
SIS with velocity dispersion a = 120 km s~^, for the group an NFW profile with c = 5 and r200 = 600 kpc. The red dashed lines 
show the excess surface densities AS of the galaxy SIS profile, the green dotted lines show that of the group NFW profile and the black 
solid lines show the sum of both. The upper left panel shows the profile for a galaxy directly in the group centre, the upper right, middle 
and lower left panels show the profile of a galaxy at distances of AD = 7, 100, 200 and 400 h~^i:pc from the group centre. As can be 
seen the role of the group contribution to AS strongly decreases with increasing distance. The lower right panel shows the mean excess 
surface mass density AS for an ensemble of galaxies distributed homogeneously within a disc of 600 kpc radius. For the mean excess 
surface mass density the contribution of the group halo is negligible out to a scale of 200 kpc. 
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scales, then on average the group acts as a mass sheet which 
does not enter the mass excess density). We conclude that 
the galaxy environment has little impact on the measured 
excess surface mass density on scales up to 200 kpc 
unless the galaxy haloes are limited to a very restricted 
area within their parent halo, i.e. if they would exist as 
central galaxies only (in this case we would measure the 
'group' excess surface density, potentially modulated by 
a miscentring term) or if they would preferentially live in 
the outskirts of group haloes. On scales larger than the 
extent of the galaxy distribution the measured combined 
signal approximately equals that of the 'parent halo'. I.e. 
the signature of satellite galaxies equals that of the isolated 
galaxies on small scales and on larger scales has a bump 
where the precise scale and amplitude depends on the 
(mean) properties of the central halo these satellites reside. 



2.5 Maximum likelihood analysis 

ISchneider fc Rixl l|l997ll first applied the maximum likeli- 
hood method to GGL to infer halo parameters. Follow- 
ing their work we will compare the measured (background) 
galaxy ellipticities to the predictions of parameterised pro- 
files. The profile parameter values are tied to that of a ref- 
erence galaxy with luminosity L* , the parameter values for 
other luminosities are determined by the scaling relations we 
have introduced in Section 12.31 The log-likelihood is given 

by 



p'Y „model 



(37) 



where aj are the two components of the polarisations 
for the j-th galaxy and, P~' is the shear polarisability, gij 
are the analytic shear values for a given profile and CTc 
quantifies the shape noise and measurement errors. The 
definition and a description of the polarisations is given in 
Section [3.41 in more detail (see equations 1411 to 149 1) . 

The advantage of the maximum likelihood analysis is 
that 'neighbouring' haloes in transverse separation (corre- 
lated structures like 'parent' haloes or satellite galaxies) 
but also along the line-of-sight (multiple deflections, see 
lBrainerd|[201 0| for a study of their relevance) are automati- 
cally accounted for correctly. 

We use the observed quantities (locations, redshifts, SED 
types and luminosities of galaxies) and the best fit halo 
parameters obtained from the maximum likelihood analyses 
to perform 3-D-projected lensing signal simulations. In this 
way we create synthetic noise-free shear catalogues from 
which we then can predict the tangential shear and the 
excess surface mass density signal and compare them to the 
signal derived from our observations. 



3 DATA 

In this Section we describe the data used for our weak lensing 
analysis. Section 13.11 gives an overview over the CFHTLS- 
Wide imaging data. In Section [3.21 we describe the spectro- 
scopic data used for the photometric redshift calibration. 



The photometric catalogue creation and the photometric 
redshift estimation are presented in Section 13.31 We give 
the details on the galaxy shape estimation in Section 13.41 
and describe the properties of our lensing sample in Sec- 
tion (33] In Section [3.61 we define our weak lensing lens and 
source subsamples. In Section [3.7l we estimate the lens mass 
errors due to photometric redshift uncertainties and shape 
measurement errors and further refer to possible biases in 
Section 13.81 due to lA. In Section 13.91 we define lens sub- 
samples with respect to galaxy luminosity, the galaxy SED 
type ('red' and 'blue' galaxy types) and galaxy environment 
density. 



3.1 Imaging data 

Our work uses publicly available CFHT Legacy SurvejQ 
Wide and Deep field data (CFHTLS- Wide and CFHTLS- 
Deep) . The surveys are descr ib ed in e.g. Hoekstra et al] 
JSemboloni et al.l (|2006D . ICoupon et al l (|2009D aiid 



11200' 



lAstier et al. (2006.) ■ They map over 190 deg^ and 4 deg^, 
respectively, in the u* , g' , r' , i' and z' broad band filters 
under superb seeing conditions of typicall y 0.8 arcsec with 
the MegaPrime camera. MegaPrime (sec Bou lade et aLll2003t ) 
is an optical multi-chip instrument with a 9 x 4 CCD array 
(2048 X 4096 pixels in each CCD; 0.186 arcsec pixel scale; 
FOV ~ 1° X 1° total field of view). The Wide and Deep sur- 
vey fields are done in 4 patches each, Wl (Dl), W2 (D2), 
W3 (D3) and W4 (D4). Coordinates and more details can 
be found at the CFHTLS- Wide webpag^fj- We restrict our- 
selves to 124 deg^ (49 in Wl, 25 in W2, 30 in W3 and 20 
in W4) of the Wide survey as the remaining ones only had 
incomplete photometric data for one or more bands at the 
time of our data reduced. We also consider the Dl and D3 
Deep survey fields as they overlap with Wl and W3, respec- 
tively, allowing spectroscopic calibration for fainter objects 
than the Wide flelds, and thus offering the possibility to es- 
timate photometric redshift accuracy for the same filter set 
as for CFHTLS- Wide but over a broader magnitude range. 
We download the preprocessed single frame data from the 
the Canadian Astronomical Data Centre (CADCfl The 
prepr ocessing with the Elixir (see iMagnier fc Cuillandrd 
[2004') pipeline includes bias and dark subtraction, flatfield- 
ing, fringe correction in the i'- and z'-band data, as well as 
photometric calibration and a preliminary astr ometric solu- 
tion. We then use the THELI -pipelin^l (see also lErben et al.l 
2005, 2009) and derive a more accurate astrometric solution, 
remap and stack the single exposures. The w ork with the 
THELI -pipeline closely follows the procedure of lErben et al.l 
(|2009l ). They had used the United States Naval Observa- 
tory (USNO) catalogues as reference for astrometric cali- 
brations. We found that an improved astrometric solution 
can be obtained when using SDSS-DR6-object catalogues 
lAdelman-McCarthv fc et al.l |2007| ) and Two Micron All- 
Sky Survey (2MASS) catalogues l|jarrett et all bood ) for 
those fields where SDSS data are not available. This im- 
provement reduces shape artefacts in the remapping to the 



2 http://www.cflit.hawaii.edu/Science/CFHTLS 

^ http: / / terapix.iap.fr / cpit / oldSite/Descart / summarycfhtlswide.html 
* http:/ /wwwl.cadc-ccda.hia-iha.nrc-cnrc.gc.ca/cadc/ 
^ http:/ /www. astro. uni-bonn.de/~mischa/theli. html 
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astrometric solution, and is relevant in the PSF anisotropy 
correction of the objects for the lensing analysis later on, 
while it is not required in photometry. We therefore gener- 
ate a stack with improved astrometric solution for the i'- 
band data used for shape measurement and keep the USNO 
based astrometric solution for photometry of all five bands. 
The data remain on one square degree tiles (with a pixel 
size of 0.186 arcsec), i.e. we do not coadd data in the small 
overlapping areas of adjacent pointings. We instead anal- 
yse the data square degree by square degree l ater on. The 
depth s of the stacked images are like those in lErben et al.l 
l|2009h (where 37 deg^ of CFHTLS-Wide were stacked), i.e. 
the limiting AB-magnitudes are about 25.3, 25.6, 24.5, 24.6 
and 23.6 {5a detection within a 2 arcsec aperture for a point 
source) in the u*,g',r',i' and z'-bands. The i'-band PSF- 
width is close to 0.8 arcsec for all frames. The proc edure for 
obtain ing bad area masks is identical to that of Er ben et al.l 
and allows us to identify the area where good pho- 
tometry and shape measurements can be obtained. After our 
data end-processing we thus have the stacked science data, 
their error frames (as weight images) and image masks. 



3.2 Spectroscopic data 

We use public spectroscopic data from the Visible Multiob- 
j ect Spectrograph (VIM DS) VLT Deep Survey (VVDS)-Deep 
jLe Ffevre et a l.' '2005'). VVDS-F22 (|Le Fevre et all |2004| . 
l2005l : iGariiiTet al. 2008) and the Deep Extragalactic Evo- 
lutionary Probe-2 (DEEP-2) programme (|Davis et al.ll2003l . 
l2007l : IVogt et al.ll2005l : IWeiner et al.ll2005l ). This provides us 
with 3 562, 7 986 and 3 746 high quality redshifts in the Wl, 
W3 and W4 fields. From these spectroscopic redshifts 1 548 
in the Wl, 3 960 in the W3 and 3 573 in the W4 describe 
objects with i' ^ 22.5. For the W2 we have also reduced 
spectroscopic VIMOS data obtained with the low resolution 
LR-Blue and LR-Red grisms from the ESO Programme ID 
082.A-0922(B). The data reduction was carried out as de- 
scribe d on the zCOSMOS release webpag^fl and iLillv et al] 
|2003). We obtain 960 high quality redshift estimates in the 
W2-field, 944 of those have i' ^ 22.5. In total we have 10 
025 spectroscopic redshifts with i' ^ 22.5 and additional 6 
229 redshifts with 22.5 < i' € 24.0. 



3.3 Catalogue creation and photometric redshift 
estimation 

As next step we convolve for every square degree tile all 
data to the PSF which is the largest among all filters. These 
are often the M*-band filter data which have a point spread 
function (PSF) FWHM between 0.63 to 1.22 arcsec, with 
a median of 0.9 arcsec. We convolve the (better) data with 
a Gaussian. Strictly speaking this procedure is not enough 
to map the data to the same PSF, because neither the 
'good' nor the 'bad' data have a truly Gaussian PSF. For 
this reason convolving the good data with a Gaussian of 
FWHMl^^ - FWHW^^^^ does not always give the best 



search for the kernel which maps the 'good' PSF to the exact 
form of the 'bad' PSF. This is a delicate procedure because 
the PSF-size a nd anisotropy also varies over the field (see 
also lDarnelfet a l. 2009 or Hildc brandt et al.ll2012h . Luckily 
the variation in PSF size over the field is similar for the dif- 
ferent filters at least for CFHTLS data. Therefore we follow 
a compromise: We test-convolve the good data with a se- 
quence of Gaussians to match bad PSF data, and measure 
the colours of stars in different apertures in a range of 8 to 
18 pixels or 1.5 to 3.3 arcsec diameter. The PSFs are consid- 
ered to be well matched in practice when the colours of stars 
are independent of the aperture size, and the convolution 
width where this is achieved best is taken for the real con- 
volution. This procedure does not make the two-dimensional 
PSF shapes of the good data identical to the bad data, but 
it ensures that the PSF is similar enough to allow matched 
aperture photometry. After the PSF homogenisation we do 
PSF photometry within an apertu re of 8 pixels (1.5 arcsec) 
diameter with SExtractor [j (see iBertin fc Arnout3 Il996l ) 
in dual-image mode, including the weight frames in the de- 
tection and extraction frame. The objects are detected in 
the unconvolved i'-band frames with a S/N threshold of 
2(7 on at least four contiguous pixels. We make use of the 
SExtractor option to convolve the data with a Gaussian 
before detection on the unconvolved i'-band. The full width 
of half maximum of this pre-detection convolution is cho- 
sen to be 0.4 arcsec and suppresses correlated noise on a 
scale smaller than the PSF. Aperture fluxes and their errors 
are extracted on the PSF matched images. The total magni- 
tudes in all filters are estimated using the difference between 
the fixed circular aperture magnitude (magapcr) and the au- 
tomatic Kron-like aperture (magauto) in the i'-band as a 
correction to the aperture magnitude in the band of interest 

magautOFiit„ ~ magapcrFiit„ + (magautoj, — magapcr;, ) . 

In these first photometric catalogues the colours of stars 
and galaxies still vary from field to field because of remain- 
ing zeropoint calibration errors and because of the galactic 
extinction. Since the CFHTLS-Wide fields are chosen to be 
off the galactic plane the extinction is rather small and does 
not vary a lot over one square degree tiles: the maximum and 
minimum extinction in all Wide fields is 0.03 and 0.14, and 
the difference between maximum and minimum extinction 
value per square degree can be up to 0.03 for high extinction 
fields and 0.01 for fields with low extinction values. We ap- 
ply one zeropoint and extinction correction value per square 
degree field by shifting the observed stella r colours to t hose 
predicted from the Pickles stellar library l|Pickleslll998l ) for 
the given photometric system. In this way we obtain cata- 
logues without remaining systematic field to field variations 
in the mean stellar and galaxy colours. 

We then derive first test photometric r edshifts. The 
photo metric redshift algorithm is that of iBender et al.l 
(|200ll ). It was successfully apphed in a variety of con- 
texts (see iGabasch et all l200 4a. 2006', l2008l: iFeulner et al l 
I2OO5I: iFeulner. Hopp fc Botzler. .2006i: jDrorv et all I2OOI 



results for PSF photometry in practice. One should instead 



^ ry 

Bri mioulle et al.ll2008l. iLerchster et al.1 [20111 . ISpineUi et all 
I20lj and recently iGruen et aLll2013l). T he SEDs used have 
been developed bv lBender et al. 1I2OOII) for h i gher redshift 
and fainter HDF (see e.g. lArnouts et al.|[l999l : ICohen et ahl 



^ http:/ /irsa.ipac.caltech.edu/data/COSMOS/spectra/z- 
cosmos/Z-COSMOS JNFO.html 



http: / /www. astromatic.net/software/sextractor 
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Figure 2. Photometric redshift accuracy (density plot) for the CFHTLS Wl to W4 fields for a depth of i' ^22.5. Each panel shows the 
object number Ng^i, the outlier rate rj, the photometric redshift scatter o"Az/(i+z) ^^'^ the mean photometric redshift error Az/{1 + z). 



I2OOOI: iFerguson. Dickinson fc Williamsll2000l') and FD F (see 
e.g. lGabasch et al.ll2004al lbl. l2006l : lFeulner et al.ll2005h galax- 
ies. Some of these SEDs actually were made to match galax- 
ies at redshifts between 3 and 4, and between 4 and 5, re- 
spectively, which are a minority in the CFHTLS- Wide data. 
We therefore take over th is SEP sample but replace a few of 
them wi th 13 SEDs fromlllbert et al. (120061 ') based on mod- 
els from lColeman. Wu fc WeedmanI ( 1980l ) to better match 
local, star-forming blue galaxies. We however do not further 



optimise the SED set for this work. We neither claim that we 
know the photometric throughput to ultimate precision (this 
is relevant mostly for the ti*-band), nor do we claim that our 
SED templates match the colours of CFHTLS galaxies per- 
fectly. We therefore use spectroscopic redshifts in our fields 
and test which empirical zeropoint offsets have to be applied 
relative to the stellar template match in order to obtain a 
good overall photometric redshift performance. For each of 
the four CFHTLS- Wide fields we determine independently 
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Figure 3. Photometric redshift accuracy (density plot) for the CFHTLS Dl- and D3-Deep and Wl- and W3-Wide fields for a depth of 
22.5 ^5 i' ^ 24.0. For these faint galaxies our photometric redshift are worse than they are for the bright sample in Fig. [2] We have shown 
in the upper panel (by using CFHTLS-Deep Dl and D3 data) how good the photometric redshifts can become for the same spectroscopic 
control sample if the photometric errors could be drastically decreased. 



the necessary offsets from the average offsets for all particu- 
lar subfields with spectroscopic data. We apply these offsets 
to the Wl to W4 fields without spectroscopic control sample. 
We have verified that this method yields reliable predictions 
by considering only a part of the spectroscopic subfields for 
the offset determination, and controlling the results with the 
remaining subfields. Also, these additional photometric off- 
sets are rather small. The offsets applied are 0.01 and 0.02 



for all filters and all four CFHTLS- Wide fields. The only 
exception is W2 and W4 where we apply an offset in the u* 
band of -0.05 and -0.06. 

We then finally obtain the photometric redshifts for 124 
deg^ of the CFHTLS- WIDE. We evaluate the photomet- 
ric redshift accuracy using all available spectroscopic data 
(see Section |3.2|I . We only consider VIMOS objects with a 
spectroscopic flag of 3, 4, 23 or 24 (confidence^ 95 per 
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cent) and DEEP2-objects with a spectroscopic flag of 4 
(confidence=100 per cent). We use our bad area flag maps 
and compare photometric redshifts and spectroscopic data 
only in areas with good photometry. Also later on photo- 
metric redshift data are only used in good area regions. 
The quality of photometric redshifts is quantified with three 
numbers: 

(i) the outlier rate rj, defined as the fraction of objects 
exceeding a rest frame error of 0.15, i.e. 

r, = fraction with ( ~ '^^^°^\ ^ q.IsI , (38) 

(ii) the photometric redshift scatter, calculated from the 
width of the central part of the error distribution 



cr^2/{n-z) = 1-48 X median or < - — ^ 



_L_ 2, I 

spec J non — outliers 

(39) 



(iii) the mean photometric redshift error 



Az/{l + z) = 



E 



1 ~h -Sgpe 



(40) 



In these equations Nspac is the spectroscopic sample size and 
Zspec and Zphot are the spectroscopic and photometric red- 
shifts. 

Fig. [2] compares photometric and spectroscopic redshifts for 
galaxies with i' < 22.5 in all four Wide fields. For fainter 
magnitudes only Wl and W3 have spectroscopic redshifts. 
Fig. [2] shows that the photometric redshift performance is 
roughly equal in all four fields. The outlier rate is between 
1 and 3.5 per cent and the root-mean-squared (rms) er- 
ror of the distribution within the dashed lines (marking a 
rest frame redshift error of 0.15) equals 0.028 to 0.034. The 
photometric redshift quality for bright galaxies will mostly 
be relevant for the lenses under consideration. Most of the 
lensed objects will be between 22.5 ^ i' ^ 24.0 however. 
For these faint galaxies spectroscopic data are only avail- 
able for the Wl and W3 fields (coming from the VVDS-Deep 
and DEEP2 data sets). We present the redshift accuracy for 
these galaxies in the two lower panels of Fig. O For these 
faint magnitudes we have outlier rates of about 5 per cent 
and an outlier-cleaned rms error of 0.04 — 0.05 (l + z). There 
are two potential reasons why the errors are increased rela- 
tive to the bright sample: first of all, the photometric errors 
are increased; secondly, at fainter magnitudes the variety 
of SED types becomes larger and the template set could 
be not flexible enough anymore or there are redshift de- 
generacies not present in bright samples because luminos- 
ity priors forbid break mismatches where the second, high 
redshift solution implies an extremely, unreasonably bright 
object. To investigate this point we also add our photomet- 
ric redshift performance for the case of the much deeper 
Dl and D3 fields. These fields have integration times of 9 
hours to 70 hours per filter compared to 20 minutes to 1.5 
hours for the CFHTLS-WIDE. We see that the faint ob- 
jects (i' = 22.5 — 24) perform better in the deep than in the 
shallow CFHTLS data sets. In fact their photometric red- 
shift scatter in the Deep fields is not larger than the redshift 
scatter of the bright, i' < 22.5 objects in the Deep fields, 



which means that our SEDs describe the bright and faint 
objects equally well, and the reduced performance for the 
22.5 < i' < 24.0-objects in the WIDE-data cannot be dom- 
inantly caused by an increased SED variety for the fainter 
objects. We conclude that the photometric error (implying 
also insignificant flux estimates in some bands) is the ma- 
jor reason why the photometric redshifts for faint objects in 
CFHTLS- Wide deteriorate. 

W e also compare the accur acy of our redshifts with those 
of lHildebrandt et all (|2012l '). who improved their CFHTLS 
data quality by executing a global homogenisation of the 
PSF over the field of view of each pointing. Fig. U shows 
the photometric redshift scatter 0-^2/(1+2), the outlier rate 
77 and the bias as a function of apparent magnitude i' and 
photometric redshift Zphot . Our results are in general similar. 
This implies that our pragmatic choice for PSF-'equalising' 
relative to a more rigorous PSF equalisation is not a limit- 
ing factor for the photometric accuracy that can be reached. 
The bias for our faintest luminosity bin is poorly estimated 
due to low number statistics. 

We classify our catalogue objects with respect to size, photo- 
metric redshift uncertainties and saturation. Saturated ob- 
jects are entirely excluded in the later analyses. Further ob- 
jects with redshift uncertainties Azphot > 0.25 (1 -I- Zphot) 
and objects with sizes smaller than the PSF are excluded 
from the background sample. 

How photometric redshift errors may bias the lens surface 
mass density estimates is discussed in Section [3.71 



3.4 Shape and shear estimates 

This section describes how we measure object shapes and 
ob tain a gravitational shear estimate. We use the method 
of Kaiser. Squires fc BroadhurstI l| 19951 . KSB), extended by 
iHoekstra et al.l 1 19981 , KSB-I-) and refer the reader t o the 



Bartelmann fc Schneider 



Bartelmann fc Schneider! (|200ll ) and 



(2001). We 



detailed derivation in 
follow the notation of 
summarise only the most important equations here. 
Object ellipticities are defined in terms of second moments 
of the surface brightness distribution I{0), 



Qii — Q2 



2iQi2 



Q11 + Q2 



(41) 



where 



Q,j = J cPe Wr^{\e\) i{e) . (42) 

Ellipticities (and gravitational shears) are complex quanti- 
ties describing the size and the major axis direction of the 
ellipticity (and size and direction of gravitational shear). The 
isotropic weight function Wrg is a Gaussian, with a width 
rg adjusted to the object size. We choose rg to be equal to 
the SExtractor flux radius in our work later on. 
The apparent ellipticities of objects are altered by the re- 
duced gravitational shear 17 = 7/(1 — k), by the PSF smear- 
ing (the isotropic part of the PSF, which makes galaxies 
larger) and by the anisotropic part of the PSF (which dis- 
torts objects). The relation of these quantities was derived 
by KSB in the linear regime of reduced gravitational shear, 
nearly Gaussian PSFs and small PSF anisotropics: 



^obs — f^s ^~ 



(43) 




Figure 4. Comparison of the photometric redshift ac curacy. The red d a shed Unes show the accuracy of our photometric redshift 
estimation, the green solid hnes show the results from iHildebrandt et aL I hOl^ ) on CFHTLS data. The plot shows the photometric 
redshift scatter crAz/(l+z) (upper panels), the outlier rate jj and the bias as a function of apparent magnitude i' (left panels) and 
photometric redshift ^phot (right panels). The results are in general quite similar. There are no objects left in the source sample with 
magnitudes fainter than j'=24 (dashed line in the left panels). The bias for our faintest luminosity bin is poorly estimated due to a few 
outliers in a very small number sample. The dashed line at redshift z=l in the right panels shows the redshift limit for our lens sample. 



In this equation eohs is the observed object ellipticity, e* 
would be the observed ellipticity of a star at the same po- 
sition, Es is the ellipticity the object would have if there 
was only a PSF-smearing (circular PSF and no gravitational 
shear), g is the reduced gravitational s hear and is the 

f re-se eing polarisability introduced by iLuppino Sz Kaiseil 
1993), 



rysm 

pg _ psh _ psh>^ 



(44) 



The second order tensors P""^ and P'*^ describe the reac- 
tion of the surface brightness distribution of an object to 
smearing and shearing by the PSF and the gravitational 
field. They are calculated from the surface brightness distri- 
butions of the galaxies and stars, where quantities measured 
for stars carry an asterisk. The division by P""^* in equa- 
tion (|44|l is strictly speaking a tensor inversion, but often 
replaced by dividing by half of the trace. In equations (|43p 
and (144 |l the values for the stellar quantities at the locations 



of galaxies have to be used. They are estimated by using 
stars in the same field of view and fitting a two-dimensional 
polynomial to describe the spatial variation of these stellar 
quantities. 

We now define the corrected ellipticities as 



^obs 



1 



(45) 



The assumption that sources are randomly oriented on the 
sky implies that their mean ellipticity and also their mean 
PSF-smeared ellipticity is zero, i.e. 



(Cs) = , 



and thus 



(ecorr) = g : 



(46) 



(47) 



holds. In the weak lensing case k ^ 1 the equation simplifies 
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to 

(ecorr) = 7 ■ (48) 

We are usin R a KSB+-implementa tion adapted from the 
TS-pipeline ijschrabback et al.ll2007l ). which was kindly pro- 
vided by Thom as Erben an d Tim Schrabback, mostly us- 
ing code from lErben et al.l (|200ll ). which itself is based 
on Nick Kaiser's original IMCAT toolfl The TS- pipehne 
has been tested in t he Sh ear TEsting Programme (STEP, 
see iHevmans et alJ l2006l ). The analysis of the first set 
of image simulations (STEPl) showed a significant bias 
l|Hevmans et all 120061 ) which could be mostly eliminated 
with a shear calibration factor c^ai = 1/0.91, 

7 = Ccal (ecorr) ■ (49) 

The analysis of STEP2 (|Massev et al.ll2007l 'l seemed to in- 
dicate that this bias calibration is accurate to ~ 3 per 
cent. The GRavitatio nal lEnsing Accur acy Testing 2008 
(GREAT08) team fsee iBridle et al.ll2010l ') showed that the 
multiplicative shear bias can be still of the order of almost 
5 per cent (at least for KSB in the Heymans implementa- 
tion), and that this bias is present for the low signal to noise 
objects, {S / N) Q-^-^^r^^^ K, 10. However, the correct S/N for 
these GREAT08 objects with {S/N)^^^^^^.^,^^^ = 10 is more 
like [S / N) Q-^^j^r^Q^_^^^^ = 6. We therefore conclude that 
even with the bias correction in equation H49p we still might 
have a bias of 5 per cent for objects with S/N « 5. The size 
of the bias can however be tested later on by comparing the 
shear signal of low and high signal to noise objects in the 
background of the same foreground structures. The same 
test can be made to analyse whether the smallest objects in 
the shape catalogue have biase s relati v e to t he larger ones. 
According to the work of Bridle ^t al.l (|2010l ) we do not ex- 
pect a bias if the ratio of galaxy and stellar FWHM exceeds 
1.4, but we expect a bias of about 7 per cent in the case 
where this ratio is only 1.2. 

To obtain the shape catalogue we run the SExtractor soft- 
ware with a signal to noise detection threshold of Scr for four 
contiguous pixels. We then measure the objects' shapes and 
their smear and shear tensors. Then bright but unsaturated 
stars with magnitudes in the range of 18 < i' < 22 are se- 
lected, where the exact limits depend on the S/N and PSF 
of each CFHTLS-field. Further we require a SExtractor star 
classification of greater than 0.96. These stars are taken to 
measure the PSF anisotropy pattern. Examples are shown 
in the upper panels of Figs. [5] and [6l where we have drawn 
a line at each location of a star, with length proportional 
to the anisotropy and the orientation along the major axis 
of the anisotropy ellipsoid (a so-called whisker plot). As can 
be seen in these Figures we do not use stars (and galaxies) 
which come from different chips in the individual exposures. 
This results in 36 distinct regions from which shape esti- 
mates are used. The gaps between these regions depend on 
the dither pattern which was originally meant to fill the chip 
gaps. Since the camera is an array of 9 x 4 CCDs, and gaps 
in vertical direction are small, one can hardly see the verti- 
cal gaps in Figs. [5] and [6l The reason why we ignore coadded 
data that comes from different chips is because there are dis- 
continuities in the PSF anisotropy fields over chip borders. 

* http://www.ifa.hawaii.edu/~'kaiser/imcat/ 



The coadded data in these regions are a superposition of 
two or more anisotropy patterns. It is impossible to obtain 
a meaningful PSF anisotropy pattern from stars in these re- 
gions and to correct the objects accordingly. Our procedure 
reduces the number of galaxies in our sample, but also the 
remaining systematic errors of the PSF anisotropy correc- 
tion. 

Ignoring these regions we then take a fifth order polyno- 
mial to describe the PSF anisotropy pattern over the one 
square degree tiles. The fifth order is flexible enough for 
the complex PSF anisotropy patterns. We have also experi- 
mented with third order polynomials for subfields of the one 
square degree fields but obtained larger systematic errors 
in terms of PSF anisotropy model residuals and B-modes 
in the two-point shear correlation fun ction (for a definition 
and measurements of B-modes see e.g. IVan Waerbeke et al.l 
2000, 2001] or lFu et al.ll2008l '). A sophisticated analysis that 
provides an objective measure whether to use lower or higher 
order pol ynomials in t he PSF anisotropy correction is de- 
scribed in lRowd l|2010l 'l. 

We show the residuals of the PSF anisotropy model and the 
data in the upper right panels of Figs. [5] and (6] In the lower 
left and right panels the ellipticity components ei and £2 of 
the stars before and after correction are shown. Their av- 
erage values and their dispersion as well as the number of 
stars used for these models (1 603 and 1 521) is given in 
these Figures as well. The data shown in Figs. [S] and [HI have 
been taken before and after a lens flip (the orientation of 
lens L3 has been inverted) at the MegaCam camera optic^. 
The lens flip significantly changed the image characteristics 
of CFHTLS observations. 

This is why the anisotropy pattern is so remarkably differ- 
ent. In general, the large anisotropy in the early data is more 
difficult to correct for than the less complex anisotropy pat- 
tern in the later data, and several frames of the pre-lens flip 
data have a much worse anisotropy residual than the exam- 
ple shown in Figs. [5] For this reason we decided to discard 
35 CFHTLS fields which performed distinctly worse after vi- 
sual inspection and only consider 89 fields for our analysis. 
After the PSF anisotropy modelling is done the objects are 
PSF anisotropy corrected. For this we calculate the stel- 
lar anisotropy correction quantities for a weight function 
with a width of rg that equals the half light radius of the 
object, Tg = rh, where the half light radius rh comes from 
SExtractor . 

We show the selection of stars (red) and galaxies (green) 
in the magnitude-half light radius plane for the field 
Wlm2m2 for illustration in Fig. [T] The resulting his- 
tograms for the FWHM of stars and galaxies can be 
seen in Fig. [SI which shows that most of galaxies obey 
FWHMgai > 1.4 X FWHMstar ~ 0.84 arcsec, which is the 
limit where no bias for the shear measurement signal was 
found in the KSB-implementa tion of Heymans according to 
GREAT08 (|Bridle et al.lbOld l. In addition when the shape 
and photometric redshift catalogue are merged to obtain the 
background galaxy sample all objects with small FWHM are 
eliminated (see dashed histogram in Fig. [8]). 
Finally, we exclude objects which have u = {S/N)-^^-^ ^ 5 



9 http://www.cfht.hawaii.edu/Science/CFHTLS- 
DATA/cfhtlsgeneralnews.html#0007 
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Figure 5. PSF anisotropy for an early (pre-lens flip) CFHTLS fleld (Wlp2p3). The upper left panel shows the observed PSF anisotropy 
pattern, the upper right panel shows the residual stellar anisotropy after modelling the two-dimensional PSF anisotropy pattern with 
a flfth order polynomial over the whole field of view. The amplitude of the anisotropy is given by showing the length of a 10 per cent 
anisotropy on the upper left of the two upper sub panels. The lower left panel shows the stellar ellipticity distribution before the PSF 
anisotropy correction, the lower right panel shows the stellar ellipticity distribution after the correction. The mean stellar ellipticity 
components, their dispersions, and the number of stars used for the modelling are also added in these lower sub panels. 



because for those shape measurements become very noisy 
and shear estimates are potentially biased low. We'll inves- 
tigate later when measuring the GGL signal whether this 
cut is justified. There is no common agreement above which 
S/N galaxies should be used for shear measurements. It cer- 
tainly depends on the science case as well. Umetsu (priv. 
communication) suggests in the context of measuring the 
weak shear profile in the outskirts of clusters to cut all ob- 



jects with S/N smaller than 7. 

We call the catalogue that remains after the KSB pipeline's 
and the S/N and size cut the 'shape catalogue'. The distri- 
butions of the two ellipticity components ei and £2 are shown 
in Figs. [51 and 1101 as green solid lines. The distributions are 
very similar and have widths of CTej = a^^ = 0.29. We also 
show the distributions for those objects in the shape cata- 
logue that also have a photometric redshift estimate as green 
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Figure 6. Same as Fig. [5] but this time using data for the field Wlmlm2 which was observed after the lens flip. The PSF anisotropy 
pattern is more regular and also has a smaller amplitude. 



dashed lines. These distributions have the identical widths 
of again 0.29; these values are used for the error analysis 
later on. The galaxy ellipticity distributions however have 
broader wings than Gaussians. In Fig. [TT] we show the dis- 
tribution of the absolute value of ellipticities as solid (shape 
catalogue) and dashed (after merging shape with photomet- 
ric redshift catalogue) curves; its medians are = 0.31. 
The distributions of the half light radii and the S/N of the 
galaxies can be seen in Figs. [12] and 1131 The median source 
half light radius is 0.6 arcsec and the median S/N-ratio for 
the objects in the shape catalogue is S/N^^d ~ 11.90. Also 
in this case the numbers remain the same when the shape 



catalogue is merged with the photometric redshift catalogue. 
This is because most of the objects in the shape catalogue 
do have a photometric redshift, whereas there are many ob- 
jects in the photometric redshift catalogue that have a too 
low signal to noise to enter the shape catalogue. 

3.5 Properties of galaxies in the photometric 
redshift and shape catalogues 

The i'-band detection based, PSF-matched total multi- 
colour aperture photometry catalogue for the 89 deg^ anal- 
ysed comprises 12 060 104 objects. The number of galax- 
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Figure 7. Distribution of the i'-magnitudc vs. SExtractor half 
light radius for objects detected with SExtractor for shape mea- 
surements in the WlmlmO-field. We mark stars in red and objects 
(galaxies) that enter our shape catalogue in green. 



ies in our photometric redshift catalogue is reduced relative 
to the full photometric catalogue: a fraction of objects is 
eliminated because not all area is appropriate for precise 
photometry (areas around bright stars, satellite tracks and 
image artefacts have to be flagged), another fraction of ob- 
jects is identified as stars by their morphology (SExtractor 
classification) or by their colours (SED match to stellar li- 
brary). 

Our photometric redshift catalogue contains 8 315 162 galax- 
ies. Their number counts are shown in Fig. [ill They become 
incomplete at i' ~ 23.5. The galaxies' photometric redshift 
distribution is shown in Fig. 1151 The majority of galaxies 
is below redshift one. The galaxy density strongly declines 
towards z=2, and has a long redshift tail up to 2 = 4. A 
few of these high redshift galaxies are also spectroscopically 
verified as can be seen in Fig. |3l O ur photometric redshift 
dist ribution look s simil ar to that of iBrimiouUe et aD (|2008l ) 
and lErben et al.l l|2009l ). who used CFHTLS-Wide data as 
well, however with a smaller field (about 37 deg'^) and a 
slightly different selection criteria for the galaxies in these 
previous works. 

Those galaxies which have photometric redshifts and shape 
measurements define our maximum background sample. 
They are 2 960 048 ahogether, where 2 416 426 reside in ar- 
eas that are not flagged as bad area. These are the galaxies 
that we take for the weak lensing analysis in the following. 
We include the number counts and photometric redshift dis- 
tributions for galaxies with photometric redshifts and shape 
estimates in Figs. [14] and [15] as dashed curves. The differ- 
ences to the (cleaned) photometric redshift sample are ex- 
plained as follows: Very bright and large and thus nearby 
objects do not end up in shape catalogues which reduces the 
galaxy density at low redshifts in shape catalogues. On the 
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Figure 8. Distribution of the FWHM of stars (red) and galaxies 
(dashed green) detected for shape measurements in the WlmlmO- 
field. The black dotted line shows the complete galaxy sample, 
while the dashed line shows the cleaned and selected galaxy sam- 
ple for shape measurements (S/N^s 5 and Zpi^^t ^ 2). Most of the 
galaxies in this field have a FWHM that is larger than 1.4 times 
the FWHM of stars (0.6 arcsec). This has been shown to allow a 
bias free shear estimate in GREAT08 ll Bridle et a"l]|2010h . After 
the photometric redshift and shape catalogue are merged there is 
another selection made which ensures that the FWHM of 'lensed' 
objects is larger than of high S/N stars in each of the CFHTLS 
frames. 



other hand, very small (close to point like) objects are elim- 
inated as well as faint objects because galaxies have to be 
larger than stars in terms of half light radius and one needs 
a higher S/N for shape measurements than for photometry. 
This removes 70 per cent of all objects in the photometric 
redshift catalogue from the shape catalogue. Faint objects 
are sp read over a broad redshift range. Fig. 14 in lllbert et al.l 
(|2006l ), for instance, shows that the redshift distribution of 
23 < i' < 23.5 and 23.5 < i' < 24 galaxies has a similar 
shape for z < 2. This implies that the shape of the redshift 
probability distribution is hardly altered by our 5'/A'' > 5 
cut in the galaxy shape catalogue relative to the photometric 
redshift catalogue although of course the number of galaxies 
is drastically reduced by this cut. For comparison we also 
have added the photometric redshift distribution for galax- 
ies in the weak lensing sample (galaxies with shape and pho- 
tometric redshift information) of ,Fu et al., (2008i). The red- 
shift distributions look rather similar in general, but due to 
stricter magnitude requirements for our source sample our 
number counts are slightly lower above z ~ 0.8. 
The SED template fltting provides an absolute magnitude 
estimate for each object at its photometric redshift. We show 
the My.' — z relation for the photometric redshift sample in 
Fig. 1161 Galaxies fainter than M^' — —18 are limited to 
z < 0.6, whereas galaxies with Af^/ < —20 can be found out 
to 2: = 1. 
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Figure 9. Distribution of the ei component of galaxies in the 
shape catalogue (solid green) and after merging the shape cata- 
logue with the photometric redshift catalogue (dashed red). The 
rms-widths in both cases are equal to ctjj = 0.29. 



Figure 11. Distribution the absolute seeing-corrected ellipticities 
|e| of objects in the shape catalogue (solid green) and after merg- 
ing the shape catalogue with the photometric redshift catalogue 
(dashed red). The median value equals 0.31 in both cases. 
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Figure 10. Distribution of the £2 component of galaxies in the 
shape catalogue (solid green) and after merging the shape cata- 
logue with the photometric redshift catalogue (dashed red). The 
rms-widths in both cases are equal to CTej = 0.29. 
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Figure 12. Distribution for galaxies in the shape catalogue (solid 
green) and the shape catalogue merged with the photometric red- 
shift catalogue (dashed red) ; the vertical line indicates the median 
of ''h mcd =0.6 arcsec. 
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Figure 13. S/N-ratio-distribution (we take the t'- value of the 
KSB-pipeUne as the S/N estimate) of objects in the shape cat- 
alogue and in the shape catalogue merged with the photomet- 
ric redshift sample. The vertical line indicates the median of 
S/A^njcd = 11.9. Note that the shape catalogue itself only includes 
objects with a S/N > 5. 
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Figure 14. Histogram of the i'-band magnitude for all 8 315 162 
objects in our cleaned photometric redshift catalogue (red solid 
line) and for those in the photometric redshift catalogue which are 
not eliminated in the shape estimation procedure (green dashed 
line), i.e. for all 2 416 426 objects which enter our weak lensing 
catalogue. Objects that are large and bright and those which have 
low signal to noise do not enter the shape catalogue. 
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Figure 15. Histogram of the cleaned photometric redshift sam- 
ple. The red solid line shows the full photometric redshift sample, 
the green dashed line shows the photometric redshift distribution 
of our shear catalogue. The gap and bump at z ~ 1.6 are a feature 
of very faint objects (i ^ 24.5) and do neither affect our back- 
ground sample (we have no background with i' ^ 24) nor our 
foreground (we only consider lenses with z ^ 1). For comparison 
we show the photometric redshift distribution of lFu et al.l l|2008l ) 
(black dotted line). 



3.6 Defining lens and background galaxy samples 

The procedure for defining a 'foreground' and 'background' 
galaxy sample or a lens and lensed object sample is guided 
by the following facts: We want to study the mass distribu- 
tion of galaxies as a function of redshift, rest frame lumi- 
nosities, SED-types and later on also stellar masses. Hence 
we can only consider galaxies for which we have photomet- 
ric redshifts. Among those we restrict the analysis to those 
galaxies from the photometric redshift sample with 



lens < — 17 and zions ^ 1 



(50) 



This defines our maximum lens sample. We can study the 
properties of any lens subsample, provided we have enough 
galaxies with shape measurements behind the lens sample 
considered. 

We construct our maximum background sample from those 
galaxies which have reliable photometric redshifts and shape 
measurements. We ignore sources with Zphot > 2 as the num- 
ber counts strongly drop for higher redshifts (see Fig. llSp . 
Another reason is that the accuracy of photometric redshifts 
deteriorates for 2 < Zphot < 3 without NIR information at 
the depth of CFHTLS-Wide and that the fractional con- 
tribution for Zphot > 3, where the photometric redshift ac- 
curacy improves due to u-band drop-outs, is negligible. We 
also ignore objects with large photometric redshift errors. In 
summary we require 



< 0.25 (1+Zph, 



ot sourcc/lcns. 



and 



< 2.0. 

(51) 
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Figure 17. Flowchart for our lens and background definition. 
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Figure 16. Density plot for the absolute luminosity (M^/) red- 
shift {z) relation for the cleaned photometric redshift catalogue. 
The cleaning regards masked area, stellar contamination and non 
precise photometric redshifts and is explained in the text. The 
number density of very bright galaxies at low redshifts is non- 
zero, but rather low, leading to an apparent gap in this region in 
the plot. 



This provides us with 4 942 433 galaxies in the maximum 
lens and 1 684 290 galaxies in the maximum background 
sample. 

In the following we always assign a background sample dy- 
namically to the foreground sample: for a fixed foreground 
sample we analyse the shapes of all galaxies within an an- 
gular or physical radius in the background of every lens 
galaxy. As minimum and maximum angular scale we choose 
5 arcsec and 15 arcmin. This corresponds to 3.3 kpc 
and 600 h'^ kpc ak. z = 0.05, 20 h'^ kpc and 3.8 h'^ Mpc 
at 2 = 0.5, 28 h~'^ kpc and 5.0 Mpc at z = 1. The 
GGL signal will only be evaluated out to a distance of 
2 Mpc later on. The outer angular cut-off of 15' is cho- 
sen for computational reasons. It equals a transverse sepa- 
ration of 2.7 Mpc at a redshift of 0.3. 
As the criterion for a galaxy to be a background galaxy we 
require that the photometric redshift estimate for the galaxy 
from the maximum background sample exceeds that of the 
lens by the - Zd ^ ^4Ai^ + 4Ai^ VSA^ . This photo- 
metric redshift requirement becomes equal to Zs ~ Zd ^0.1 
for our photometric redshift errors of Az ~ 0.04. 
On the other hand, even if the photometric redshifts were 
extremely precise one would not consider galaxies in the ulti- 
mate background of a lens as they only carry a small lensing 
signal. The gravitational shear observed in the background 
of a singular isothermal sphere scales as Dds/fs, where Dds 
and Ds are the diameter distances from the lens to the source 
and from the observer to the source. This ratio approaches 
zero for Zs ~ z^, rises steeply for Zs > z^ and then flattens for 
larger redshifts. To exclude galaxies that carry no signal but 



only shape noise, for any given foreground sample we reject 
galaxies for which Dds/Dg > 0.1 for a given foreground sam- 
ple does not hold for any of the galaxies in their foreground 
(cf. Fig, lisp. Finally, when the measured shear is translated 
into a mass density, the fractional error is proportional to the 
ratio of the error of the critical surface mass density and the 
critical surface mass density (see equal ion [59]). For Gaussian 
photometric redshifts errors with a width of 0.05 {1 + z) this 
ratio is smaller than 0.3 if > 1.1 Zd -1-0.15 and Zd > 0.05 
(see Fig. I19|l . This mass accuracy condition is slightly more 
conservative than z^ > Zd + 0.1. Furthermore, this criterium 
implies D^JD^ ^ 0.75 for Zd ~ 0.05 and DdJD^ > 0.15 for 
Zd ~ 1.0 which is the largest lens galaxy redshift that we 
consider in this work. In summary, our selection criterium 
for the lens and source redshifts are (see also Fig. I17|l : 

0.05 < Zd^ 1.0, Zs ^ 2.0 and z^ > 1.1 Zd + 0.15 . (52) 

These criteriae insure that foreground and background 
galaxies can be disentangled, that Dds/Ds > 0.1 holds and 
that the fractional statistical error for the critical surface 
density due to photometric redshift errors is smaller than 
0.3 per foreground-background pair (see Fig. [19}. 
The definition of our maximum lens and background sam- 
ples implies that the maximum lens sample is larger than the 
maximum background sample. It allows situations where the 
apparent magnitude of a foreground galaxy is larger than 
the magnitude of a background object tracing its gravita- 
ti onal shear. Th i s is n ot the case in GGL analyses like that 
of jParker et al.l (|2007| ). and de facto ^although not by de f- 
inition) not the cas e for the work of Fischer et aU (|2000t ). 
iMcKav et aU (|200ll ) or|Mandclbaum ct al. (200^U, due to 
the extremely bright foregroun d sample (r' < 17.7 7 for the 
main spectroscopic sample, see IStrauss et al.|l2002l and r' < 
19.5 for the spectroscopic LRG sample, see Eisenstein et al.l 
I2OOII ). Also, in our analysis the redshift distribution of lenses 
and lensed objects overlap, although for each lens only galax- 
ies in its background are cons i dered. 

Similar to lMandelbaum et al.l (|2006d ) we split the maximum 
lens sample in absolute magnitude intervals of one magni- 
tude width from M^, = -17 to M^, < -24. 
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3.7 Lens mass errors from photometric redshifts 
and shape measurement errors 

As an estimator for equation ([6]) we use tlie foreground- 
background pair average, i.e. 



AE(J?) = (Ee 7t(^) 



fg — bg — pairs 



(53) 



In the upper equation the estimates carry a 'hat'. Ec is 
obtained from the photometric redshift estimates of the 
background- foreground pair and the shear estimate 7t(-R) 
is estimated from the background shapes and the lens pho- 
tometric redshift because the observed angular scale 9 needs 
to be translated into a length scale with R — ODd- The re- 
lation between the true gravitational shear 7t(-R) at radius 
R and its estimate ^t{R) at radius R reads as 



7t (ii) = 7t (^) + [7t {R)-lt{R)\+ H 



shape 



; 7t(i?) + A7t(i?, R) + 57shapo 



(54) 



In equation (|54p we have introduced a radius independent 
shear estimation error 57shapc that comes from shape mea- 
surement error, intrinsic shape noise and potential shape 
estimate bias. 

The other quantity which we have introduced in equa- 
tion H54I) is the 'profile error' 

A7t(i?,i?)=7t(i?)-7t(i?)=7t(e-Dd)-7t(e-Dd) , (55) 

which comes from mixing physical scales when translating 
angles into length scales. The latter is zero when a spec- 
troscopic lens sample is used. In general, the profile er- 
ror depends on the profile steepness. For a power law with 
7t(-R) oc R~" it equals 



A7t(i?,i?) =7t(i?) 



7?" 



(56) 



Thus the profile error is the more important the steeper 
the shear profile is. For an isothermal profile with a = 1 
equation (|56l) becomes 



A7fS(7?, R) 



7t(ii) 



Dd 



For small redshifts the profile error equals 



Jt{R) 



Zd ~ Zd 
Zd 



(57) 



It is obvious that photometric redshift errors without a bias 
can already impose a bias on the shear profile estimate. 
To understand this we consider foreground-background pairs 
with a transverse distance estimate R in the lens plane and 
a redshift estimate of id. If the redshift estimate is too low, 
i.e. id < Zd and thus R < R, sources are scattered from 
larger to smaller transverse separations. This leads to an 
underestimate of the gravitational shear signal. Analogously 
if the true redshift is lower, the gravitational shear signal 
is overestimated. Therefore a possible bias depends on the 
redshift (and thus projected separation) distribution of the 
analysed lens sample. For a flat redshift distribution about 
the same number of galaxies are scattered down from higher 
and scattered up from lower redshifts into the bin of in- 
terest. In this case the mean projected distance of galaxies 



within a ring of a diameter of R equals the true value R. 
If the lens number increases as a function of redshift more 
galaxies with Zd > id get scattered down to a lower redshift 
than galaxies with Zd < id get scattered up. This can lead 
to (id) < (zd) and thus (R) < (R). We performed a lensing 
signal simulation, scattering our lens redshifts by adding a 
Gaussian redshift error distribution of 0.03 (l-fz), to esti- 
mate the maximal bias in the measurement of the velocity 
dispersion due to lens redshift errors. Even in the most ex- 
treme scenario (lowest redshift lenses, i.e. steep rise in red- 
shift counts and asymmetric redshift scattering, as there are 
no lower redshift galaxies that can be scattered up) the bias 
in a is below 4 per cent. For higher redshift lenses the bias 
rapidly decreases. As the number of low redshift lenses is 
low compared to the total lens number we make no attempt 
to correct for this bias. 

For a given background-foreground distance distribution the 
bias is the larger, the larger the photometric redshift scatter 
is. For flatter proflles the bias is smaller (see equation [56]) . 
The ratio of the estimator and the true contribution for each 
pair in equation (|53[) is in linear order: 

Sc UR) ^ 1 , ^ , A-ft{R,R) ^7shapc , . 

Ec7t(i?) Ec 7t(i?) ^ 7t(i?) ■ ^ ' 

We now deflne A"E = AE + (5AE, and Ec = Ec -I- (5Ec, insert 
equation (I59p into equation (I53p and obtain to linear order 
for the error of the estimator: 



5AE(J?) = (7.(i?)Ec[i^-,-^^ , 
7t(i?) <5Ec + A7t(i?, i?)Ec + <57shapeEc 



fg-bg-pair 
bg — pair 



(60) 

Equation H60|) can be used to obtain the error estimates in 
presence of scatter and biases for photometric redshifts and 
shape estimates. Equations (|56l) . (|57|) and H58p explain to- 
gether with equation H60p why lens redshift errors are more 
severe in GGL than source redshift errors, in particular if the 
lens is very close. In this case, for Azd ~ 0.03 and Zd ~ 0.1 
this fractional error can easily approach 30 per cent. For 
larger redshifts the dependency of Dd on redshift weakens 
on, so does the associated profile error. The term (57shapo 
includes the intrinsic shape noise, the shape measurement 
error and potential systematics. The value of 57shapc/7t is 
in order of 0.3/0.002 « 150 per foreground background pair. 
So, the relative shape error exceeds the relative proflle error 
by more than a factor of 500. 

We now calculate the errors due to shape noise and photo- 
metric redshift errors more accurately. 

Regarding the shape measurement noise one can use the esti- 
mate of the tangential shear relative to a random foreground 
sample to estimate the statistical error of 7t(-R) and AE(ii) 
estimates. This is shown in more detail in Appendix lAll For 
a given lens sample as defined in Section r3.6l the correspond- 
ing background sample is specified by the selection criterion 
defined in equation (|52p . We do the same estimate of 7(-R) 
and AE(ii) as we do for real data but use the shape of the 
background object with a randomised phase or the shape 
of another object in the shape catalogue instead of the true 
tangential ellipticity. 

In addition we simulate photometric redshift errors by 
adding a Gaussian distributed error of 0.05 {1 + z). We cal- 
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culate the lens distance and the critical surface mass 
density AEc and estimate the propagated errors for both 
quantities due to the photometric redshift uncertainties. The 
simulation is designed to rather overestimate the estimated 
errors than to underestimate them. This is especially true 
for very low lens redshifts, where the scattering of photo- 
metric redshifts can lead to negative lens redshifts, which 
are not realised in observational data. 

The propagated errors comprise a systematic and a statisti- 
cal part and read as 



5Ec = A/5E, 



+ 



ST., 



(61) 



In order to disentangle their contributions we consider two 
samples of different sizes (100 000 and 1 000 objects) and 
solve the resulting system of equations (see also lGruen et al.l 
I2OIOI ). In both cases the systematic error is significantly 
smaller than the corresponding statistical error. For a red- 
shift of 2: = 0.05 the fractional systematic error in SD^/Dd 
is lower than 10 per cent, dropping below 5 per cent for 
higher redshifts (see Fig. llSp . The fractional error of the 
critical surface density (SEc/Ec is larger than for Dd due to 
multiple dependencies on the redshift (Dd, Da and there- 
fore Dds). However by applying the selection criterion from 
equation (|52[) we ensure that the systematic errors in unfor- 
tunate cases are at most slightly larger than 10 per cent (see 
Fig.[l9|. 

We now estimating the systematic errors for extreme cases. 
For a lens at z ~ 0.7 we expect that about 25 per cent of the 
background galaxies are close enough to the lens to have an 
average bias of about ~ 20 per cent and that the Ec estimate 
is nearly unbiased for the remaining 75 per cent. This gives 
an average systematic error of 0.75 x 1+0.25 x 1.2 — 1.05, i.e. 
an overestimation of (5Ec/Ec ~ 5 per cent. For the other ex- 
treme at low redshifts the fraction of lens-source pairs with 
a bias of ~ 20 per cent is only 5 per cent. This leads to a 
fractional systematic error of 0.95 x 1 -I- 0.05 x 1.2 = 1.01, 
i.e. 5Ec/Ec = 1 per cent. For the whole sample we expect a 
systematic bias of not more than 0.88 -I- 0.12 x 1.2 — 1.02, 
i.e. a maximal bias of 5Ec/Ec = 2 per cent due to pho- 
tometric redshift inaccuracies. We therefore conclude that 
the expected systematic errors due to photometric redshift 
uncertainties are small enough to be neglected. 



3.8 Systematic errors from lA 

Galaxies which are in the same structure and thus phys- 
ically connected are not randomly dis tributed in orienta - 
tion but rather intrinsically aligned jHirata et al.l l2004h . 
for instance satellite galaxies tend to be radially aligned 
relative to their central galaxies. This is why intrin- 
sic alignment (lA) is a major issue in the interpreta- 
tion o f c osmic shear data (s ee e.g. iMandelbaum et al.l 
l2006bl or iBridle fc Kin3 120071 '). The observed two-point 
correlation function for the ellipticity of galaxy pairs is 
{e,e*) = (7,7*) + (ef ef *) + (ri'^j*) + (ef-y*), where e''^ and 
e are the unlensed and lensed ellipticities and 7^ is the cosmic 
shear at redshift Zi (see e.g. IJoac himi fc Schneider 2008). 
The first term on the right hand side is the desired cos- 
mic shear signal, the second term (called II) describes the 
intrinsic alignment of two galaxies. Unless the two galax- 
ies are physically associated (i.e. they are required to be at 




"■lens 



Figure 18. Fractional systematic and statistical error SD^/D^ 
of the angular distance of the lens in presence of photometric 
redshift errors. The dashed line shows the statistical and the solid 
line the systematic error for a Gaussian redshift error distribution 
with a scatter of 0.05 (I + 2:). The systematic error in is below 
10 per cent for 2 ^ 0.05 and well below 5 per cent for z ^ 0.1. 



same redshift) this term is zero. The third term describes 
(for Zi ^ Zj) the correlation of a foreground gravitational 
shear with the intrinsic ellipticity of a background galaxy 
and is zero. The fourth term (called GI) describes the corre- 
lation between the intrinsic ellipticity of a foreground galaxy 
and the gravitational shear acting on a background galaxy. 
In GGL, however, one measures the tangential alignment, 
i.e. the cross-correlation of a background galaxy shape and 
the foreground lens position. Therefore intrinsic alignment 
theoretically should not be an issue at all. This situation 
is different in case of a foreground-background mismatch 
due to photometric redshift errors, where the photometric 
redshift of the assumed background object is overestimated 
and the galaxy actually is embedded in the foreground struc- 
ture. If the falsely assumed background galaxy is randomly 
oriented relative to the foreground galaxy considered, then 
the shear signal is just diluted and our error considerations 
from Section 13.71 apply. If however the background galaxy 
has a preferred direction to the foreground an additional 
source of systematic error arises. If these false 'background' 
galaxies are fainter than the foreground galaxies, they will 
likely be their satellites (if associated to the foreground 
structure) and thus will on average be radially aligned (see 
lAgustsson fc BrainerdllioOSh . This then leads to a false de- 
tection of the GTsignal. 

The separation of GGL and int rinsic alignment is investi- 
gated in detail by iBlazek et al. To isolate lA from 
the lensing signal, they exploited the fact, that the contam- 
ination of the background galaxy sample with foreground 
galaxies should decrease if a more distant background slice 
is considered. They measured the excess surface mass den- 
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Figure 19. Fractional error of the critical surface density, 
5Sc/Sc, in the presence of photometric redshift errors as a func- 
tion of the lens and source redshift. The upper panel shows the 
statistical error per foreground-background pair, the lower panel 
the systematic error. The photometric redshift errors in this sim- 
ulation are Gaussian with a scatter of 0.05 (1 -I- z). The magenta 
and the red contours show (5Sc/Sc = 0.5 and (5Ec/Sc = 0.3 levels, 
respectively. Blue and green contours show the (5Sc/Sc = 0.1 and 
5Ec/Sc = 0.05 levels. The fractional systematic error is below 0.3 
if Zsource = 1-1 ziena + 0.15 (dashed black line) and zicns > 0.05, 
which is the source-lens redshift requirement set in equation I I52I I . 
The errors are in general rather overestimated than underesti- 
mated. This is especially true for very low redshift sources, where 
the scattering during randomisation of the lens redshift can lead 
to negative lens redshifts, which are not realised in the observa- 
tional data. 



sity AS associated with SDSS-LRGs using two source sub- 
samples in two redshift slices behind the lens. 
They conclude that the size of lA for their lens sample is 
small (and consistent with zero, see their fig. 3) . Their fig. 2 
shows that for all scales larger than 100 kpc the signal 
extracted for two redshift subsets agrees which implies that 
the imprint of lA on AE can be neglected. 
In our case we can infer the potential error due to lA from 
the upper left panel of Fig. IA5I in Appendix IA4I The ma- 
genta and green points show the AE values obtained for 
foreground lenses with 0.05 < z ^ 0.5 using the shear signal 
from galaxies in the redshift slices of 0.6 ^ z ^ 0.73 and 
1.01 z ^ 2. Since the contamination of the 2 = 1 — 2 
sample should be zero, the difference between the green and 
magenta points quantifies the maximal error due to lA in 
the low z background sample with 0.6 ^ z ^ 0.73. All val- 
ues agree within 2g. We therefore conclude that systematic 
errors due to lA are small enough to be neglected. 

3.9 Defining lens subsamples 

As we do not only want to analyse the properties of 
the mean lenses as a function of luminosities, but also 
want to investigate the differences based on SED and 
environment we define several l ens su bsamples. We follow 
the approach of iDahlen et al.l (|2005l ) who used the rest 
frame [B — 1/)-colours to classify red and blue galaxy 
types. We therefore estimate the {B — l/)-colours for our 
objects from the best-fitting templates in the photometric 
redshift estimation in the AB-system and define our red 
galaxy sample as all galaxies with {B ~ V) > 0.7 and 
consequently our blue galaxy sample as all galaxies with 
{B — V) ^5 0.7. We test this selection by measuring the 
(rest frame) {B — y)-colours of the best-fitting SED (at 
the spectroscopic redshift) f or the SPSS Lu minous Red 
Galaxy (LRG) sample (see lEisenstein et al ] [2001). The 
{B - V)-colour histogram of the best-fitting SDSS-LRG 
templates (Fig. I20|l shows that this assumption excellently 
holds for the LRG sample. 

We further verify our galaxy classifications by plotting 
the apparent galaxy colours and rest frame luminosities. 
We use the absolute (<?' — r')-colour vs. the rest frame 
r'-band luminosity t o discriminate red and blue galaxies (cf. 
iLovedav et"aLll2012l ). As we can see in Fig. [21] red and blue 
galaxy populate distinct regions in the magnitude-colour 
space, especially considering higher redshifts, indicating a 
low contamination rate of our re d and blue lens samples. In 
contrast to Loved av et al.l ()2012l ) we do not have individual 
rest frame colours for our lens galaxies. We thus consider 
the apparent colours vs. the absolute magnitudes. However, 
this is at least for very low redshifts a valid approximation. 
Further we plot t he {g' — r')- vs. (r' — i')-colours (cf. 
iToieiro et al.ll2013 for their red galaxy sample at moderate 
redshifts between 0.5 and 0.7). We see in Fig. [221 that 
for redshifts z ^ 0.4 our red and blue galaxy samples 
are clearly separated and also for lower redshifts our red 
galaxy sample shows few overlap in this colour-colour-plane. 

In addition we measure a relative local galaxy density 
for all galaxies in our lensing sample. For this we consider 
all galaxies within an interval of Az = ±0.2 in front 
and behind the investigated lens galaxies. We measure 
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Figure 21. Density contours for the colour-absolute magnitude-relation of our rod and blue galaxies. The black solid lines show the 
distribution of all galaxies, the red dashed lines the distribution of red galaxies and the green dotted lines the distribution of blue 
galaxies. As we can see red and blue galaxies populate distinct regions in the colour-magnitude space, especially for higher redshifts. For 
low redshifts, there is a small overlap of red and blue galaxies in the colour-magnitude-plane. 
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Figure 22. Density contours for the apparent colour-colour-relation. The black solid lines show the distribution of all galaxies, the red 
dashed lines the distribution of red galaxies and the green dotted lines the distribution of blue galaxies. In colour-space our de fined 
red and blue galaxy samples populate distinct regions when considering redshift z ^ 0.4, in a agreement with lToieiro et ah I ll2012h . For 
redshifts lower than z 0.4 the {g' — r') vs. (r' — j')-colours of red and blue galaxies overlap. 
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Figure 20. { B — V)-rest frame colou r histogram for the SDSS 
LRG sample llEisenstein et al.l I2OOII) in the AB-system. The 
{B — V)-colours are taken from our SED-templates that best 
match the SDSS photometry at the spectroscopic redshift of the 
LRG. For almost all objects the assumption (B — V) > 0.7 holds, 
just ifying the chosen ga laxy classification for red and blue galaxies 
(see iDahlen et al]|2005^ . 



the galaxy density (galaxies per area) within an angular 
distance of 30 arcsec around a lens and on the entire field 
in the previously defined redshift interval and divide the 
quantities to derive the relative local galaxy density. We 
then define following environmental lens samples: 

(i) a very dense environment lens sample (10 per cent of 
galaxies populating the densest environments), 

(ii) a dense environment lens sample (50 per cent of galax- 
ies populating denser environments), 

(iii) a low density environment lens sample (50 per cent 
of galaxies populating lower density environments) 

(iv) and a very low density environment lens sample (10 
per cent of galaxies populating the lowest density environ- 
ments) . 
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Figure 23. Flow chart for our lensing analysis. 



haloes by explicitly modelling nearby haloes traced by galax- 
ies and analyse the AE(_R)-signal on larger scales using a 
maximum likelihood analysis that parameterises the scaling 
relations as well. We confirm the scaling relations obtained 
in the first step. Then we fix the luminosity scaling relations 
and obtain the halo parameters for the 'reference galaxy' 
as a function of SED-type and galaxy environment density. 
Finally we use the obtained values for the halo parameters 
and thus the masses to estimate the remaining scaling re- 
lations for the halo parameters, i.e. the scaling of the BBS 
truncation radius s and the NFW concentration parameter 
c as a function of luminosity. The structure in our lensing 
analysis is also illustrated in Fig. 1231 



4 WEAK LENSING ANALYSIS 

The structure of this section is as follows. First we measure 
the tangential shear signal "it{R) as a function of luminosity, 
galaxy SED and environment density to obtain a qualitative 
overview about the lensing signal strength. Using photomet- 
ric redshifts we calculate the excess surface mass density 
AE(_R). We finally use the results for different luminosity 
bins to derive the scaling relations for various halo profile 
parameters (e.g. velocity dispersion o or virial radius r2oo) 
as a function of the galaxies' absolute luminosity. In a first 
step we only consider AE(J?) on small scales and treat each 
galaxy halo as being isolated and derive the luminosity scal- 
ing relations. In the second step we account for associated 



4.1 Tangential shear 

Measuring the gravitational shear as a function of angular 
scale with respect to the foreground galaxies is the most di- 
rect measurement of GGL. The use of photometric redshifts 
is in principle not required. Instead one can define magni- 
tude intervals for foreground and background galajcies and 
then interpret the shear measurements using an assump- 
t ion on the foreground and background redshift distributi on 
fsee'Hoekstra ct al."200^: lHoekstra. Yee fc Gladdersll2004l or 
[Parkor ct al. 2007). However, if gala:xies are spread over a 
broad redshift range, then physical scales are mapped into 
various angular scales and the observed signal (as a func- 
tion of angle) is a mix of various physical scales. We therefore 
make use of the photometric redshift estimates for the lenses 
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Figure 24. Tangential shear signal (filled circles) for different lens samples. The cross shear (systematic error) is shown with empty 
squares. It is consistent with zero at all radii and in all samples. The green dashed lines in the upper left panel show the 1 — cr-level of 
the remaining systematic error. The weighting scheme is given in the text. Obviously, the fainter the subsample, the less distant it is 
on average. We use the weighted mean redshift of the galaxies to fit the amplitude of an isothermal sphere on scales R 200 kpc. 
The corresponding effective velocity dispersion of the sample is then derived as well. The shear of the 'effective' isothermal sphere is 
added. The measured gravitational shear signal is isothermal for R ^ 200 kpc. The upper left panel shows the tangential shear 
signal in the luminosity bin —24 ^ M,.i < —20 for red galaxy lenses as filled red triangles (dashed fit-line), blue galaxy lenses as filled 
blue squares (dotted fit-line) and for all galaxy lenses as filled black circles (solid fit-line). The deviations from a single halo signal due 
to the neighbouring galaxy haloes on scales R > 200 kpc are most prominent for the red galaxy lens sample as this galaxy type 
is mainly populating higher density regions such as galaxy clusters. The gravitational shear caused by the red galaxy lens sample is 
significantly higher than for blue galaxies. The upper right panel shows the combined lens sample signal separated into four luminosity 
bins of —24 M^/ < —23 in magenta (crosses, dashed-dotted fit-line), —23 ^ M^i < —22 in red (filled triangles, dashed fit-line), 
—22 M^i < —21 in blue (filled squares, dotted fit-line) and —21 :C A/^/ < —20 in green (filled circles, solid fit-line). As expected, the 
signal strength and so the effective velocity dispersion of the sample decreases with decreasing brightness of the lenses. The fractional 
decrease of the signal for the fainter samples is of the expected amplitude if TuUy-Fisher or Faber-Jackson luminosity-velocity scaling 
relations are assumed (more details in the text). The lower panels show the red galaxy lens sample (left) and the blue galaxy sample 
(right) split into the same luminosity bins with the same colour scheme as for the combined sample before. Obviously the observed 
lensing signal and therefore galaxy mass is higher in every single luminosity bin for red galaxies than for blue galaxies. 
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and measure the mean gravitational shear as a function of 
the projected distance to the foreground lenses at their red- 
shifts. The quantitative interpretation is not the major goal 
of this analysis, but to see qualitatively to which distance 
the signal is dominated by the halo the galaxy is centred 
in, and when the contribution of nearby galaxy haloes be- 
comes visible. As lenses we either consider all galaxies (in 
our photometric redshift sample) or we consider subsam- 
ples in absolute r'-band magnitudes. For each radius R the 
mean gravitational shear is obtained by a weighted mean of 
shear estimates for all foreground- background pairs at this 
radius. The w eighting factors for the shear estim ates are as 
introduced bv lHoekstra. Franx fc KuiikenI (|200Cll l 



(62) 



where CTe stands for the scatter of the intrinsic ellipticities 
and (Ae^)^''^ for the uncertainty in the ellipticity mea- 
surement. Applying these weight factors we are able to 
reduce the contribution of background galaxies with large 
uncertainties in the ellipticity measurement and to enhance 
the signal-to-noise ratio. 

Weak lensing measurements also provide a simple and 
comfortable possibility to identify systematic errors, the 
cross-shear, also called B-modes. Being a conservative force 
and therefore not producing curls, the tangential shear 
should vanish if all background objects are rotated by 45 
degrees. Vanishing B-modes do not absolutely guarantee the 
absence of systematics, but their presence is an indicator 
for remaining systematic effects. 

In Fig. [24] the gravitational shear estimate is shown 
as a function of the projected distance R to the lens when 
we average over all galaxies and over all red and blue 
galax;ies, respectively. Singular isothermal shear fits are 
added as dashed curves, where the fit has been obtained by 
using separations smaller than R = 200 kpc only. For 
the distance ratio Dds/Ds required for translating the shear 
into a velocity dispersion we use the foreground-background 
pair weighted average of the individual D^s/Ds ratios for 
each foreground-background pair. The fit shows that the 
signal approximately follows an isothermal profile out to 
R = 200 kpc. For scales larger than R = 200 kpc 
the deviations are stronger for the red galaxies than for 
the overall sample, which can be explained by red galaxies 
being more strongly correlated with each other, and by the 
fact that the fraction of red galaxies residing in denser dark 
matter environments (groups and clusters of galaxies) is 
higher than for non-red galaxies. The observed B-modes 
are consistent with zero. In the other three panels of Fig. [21] 
the signals are evaluated for four absolute magnitude 
intervals of M^/ between -24 and -20 with interval widths 
of one magnitude. As expected we see a clear sequence 
from the brightes t lens sample to the fa intest. Assuming 
a Faber-Jackson (iFaber fc JacksonI ll976l l or TuUy-Fisher 
||Tu11v fc Fished Il97 7l) relation, respectively, we observe a 
picture consistent with the decrease of velocity dispersion 
for fainter lens samples. The expectation that red galaxies 
are more massive than blue galaxies for given r'-luminosity 
is confirmed by the investigation of the individual luminosity 
bins. The gravitational shear for red galaxies is significantly 
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Figure 25. Tangential slicar signal for lens samples with lu- 
minosities —24 M^i < —17 in different environments. The 
SIS-fits are obtained within a separation oi R = 200 kpc. 
The definitions of the environments are given in Section 13.91 The 
black circles (solid fit-line) show the complete lensing sample (see 
also Fig. 1241 upper left panel), blue (filled squares, dashed-dotted 
fit-line) and green (diamonds, dashed-dotted fit-line with longer 
dashes) samples show lenses in low density environments. It be- 
comes apparent that in dense environments (red filled triangles, 
dashed fit-line) and especially for very dense environment (ma- 
genta crosses, dotted fit-line) the observed gravitational shear sig- 
nal cannot be explained with only one lens. 



higher than for their blue counterparts. The gravitational 
shear of the combined galaxy sample is between that of the 
red and blue sample for every luminosity bin as expected. 
The B-modes are consistent with zero for all luminosity 
bins. 

We now also analyse the tangential shear as a function of the 
environment density defined in Section 13.91 Fig. [5S] shows 
the gravitational shear around galaxies with luminosities of 
—24 < Mr' < —17 for the four different environment classes 
defined by us. The denser the environment, the higher the 
signal. In principle the higher signal could be due to a higher 
mean luminosity in the considered luminosity interval. But 
this is not only reason as can be seen if the behavior as 
a function of distance to the galaxies is analysed: The 
gravitational shear becomes approximately equal at the 
smallest scale of i? = 30 — 40 kpc. This implies that, 
on average, the galaxies in the four environment samples 
have roughly the same central density profiles but differ in 
their outer halo profiles and experience a different impact of 
nearby haloes. As expected, the gravitational shear is lowest 
for galaxies in the lowest environment density: it rapidly 
drops down on very sort scales and becomes negative for 
larger scales, which means that the average convergence at 
radius R is higher than the averaged convergence within 
radius R, or that the larger scale environment is denser than 
the density around the galaxy within R « 100 — 200 
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kpc itself. For the low density environment (blue points in 
Fig. I25p the shear signal follows an isothermal profile (blue 
dashed curve) out to i? = 200 kpc, and approaches 
zero on larger scales. There is hardly any impact of nearby 
haloes visible in the signal (as expected if the environment 
is poor). For the high density sample the signal decreases 
similar to an isothermal sphere out to _R = 200 kpc and 
then stays constant. For the highest environment density 
sample this happens already at scales of approximately 
R — 100 kpc and the shear stays high at a level of 
0.003. At scales of 400 h'^ kpc this signal of 0.003 is the 
same as that of a group with a dark matter halo with 
velocity dispersion of 400 km s~^ if the typical distances of 
sources and lenses are used. The shape of the signal is, due 
its 'flatness', not in agreement with being caused by a single 
group halo, where the majority of the considered galaxies 
are central. This flat behavior is also confirmed in our 
3D-LOS-projected lensing signal simulations (see Fig. IB6I 
in Appendix [B|, where we see that this flatness originates 
in the multiple gravitational deflections on brighter nearby 
galaxies in the close environment. 

4.2 Measurements of the excess surface mass 
density AE(7?) 

The measured gravitational shears and the photometric red- 
shifts of the lenses and sources are now combined to obtain 
the excess surface mass density AE(i?) (see equation [S]) ac- 
cording to the estimator in equation (|53l) . further using the 
weighting factor from equation (|62[) . We fit the excess sur- 
face density with a power law. 



AE = A[R/1 Mpc]-°' 



(63) 



where an exponent of a = 1 corresponds to an isothermal 
profile. The fits include scales out to a projected distance of 
1 Mpc. A power law fit is a too simple model to prop- 
erly describe the excess surface mass density, however, the 
profile steepness of a power law is a measure of the relative 
importance of the halo hosting the galaxy itself (slope close 
to -1), and the haloes in which the galaxy host halo resides 
(flatter slope). We split our galaxy sample samples into sev- 
eral absolute luminosity intervals in M^i and analyse them 
separately (see Fig. 1261) . We find that the slope increases 
from a ~ 0.3 for luminosities — 17 > Af^/ > —18 to values of 
a ~ 0.9 for galaxies with —23 > M^' > —24. The amplitude 
of the excess surface mass density increases with luminosity 
as seen for the tangential shear signal before (Fig. I24|) . For 
the fainter luminosity bins the signal of the total galaxy 
sample is dominated by blue lens galaxies. The excess 
surface mass density for faint red galaxies is for all scales 
significantly higher than for blue galaxies, suggesting that 
the red galaxies in this luminosity range are significantly 
more massive than their blue counterparts and that faint 
red galaxies are more likely in denser environments, where 
in addition larger scale haloes from galaxies, groups and 
clusters contribute. With increasing luminosity the signal 
amplitude increases and the difference between red and blue 
lens galaxies decreases. For the brightest luminosity interval 
the excess surface mass density of the combined sample 
is dominated by the bright red galaxies as in this lumi- 
nosity range red galaxies outnumber their blue counterparts. 



We com pare our results regarding AE to lMandelbaum et al.l 
(|2006ch . who analysed the excess surface mass density for 
luminosities —17 ^ Afr',sDSS,AB > —22.5 for red and 
blue galaxy samples using SDSS data (spectroscopic 
redshifts for the lens sample, photometric redshifts for 
bright source galaxies with ^sngg < 21 and statistical 
redshift distributions for fainter sources). In contrast to 
iMandelbaum et al.l l|2006d ) our rest frame magnitudes are 
not given in the AB but in the vega system and are not cal- 
culated for a Hubble parameter of Ho — 100 km s^^ Mpc~^ 
but for Ho = 72 km s~^ Mpc~^. To adjust for these 
differences we need to apply a magnitude offset of roughly 
A mag = —0.55 to our rest f rame magnitudes. Compared 
to IMandelbaum et al.l l|2006ct ) we use a much smaller area 
but significantly deeper data. We choose the corresponding 
luminosity bins and investigate the individual results for 
AE. The signals for the blue galaxy lens sample are in good 
agreement, yet fairly noisy. 

This is not the case for red galaxies. While for the faintest 
luminosity bins we obtain higher amplitudes up to lumi- 
nosities of Mr' ~ —21 we measure lower amplitudes for 

brighter galaxies. 

Howe ver, in th eir later work Mandelbaum. Seliak fc Hiratal 
|2008| (see also iDutton et al.TioiOl ) do find a higher sig- 
nal for faint galaxies, which agrees with our results. 
A ver y good illustration of the Mandelbaum et al] 
(l2006d). iMandelbaum. Seliak fc Hiratal ll2008f) aiid 
ISchulz. Mandelbaum fc Padmanabhan (|2010l ) results 
can be found in Dutton et al. l|2010l) . fig. 1. It shows 



that the mass-to-light (or mass-to-stellar mass) ratio of red 
galaxies increases relatively moderately with stellar mass (or 
light) i n contrast to the earlier results of lMandelbaum et al.l 
(j2006g) (which ha d still relatively la rge uncertainties). 
The compilation of iDutton et all (|2010| ). Fig.l, shows that 
the mass-to-light or mass-to-stellar mass ratio increases 
rather slowly for ste llar masses with log MBtar=10.5 to 11.3 
accord ing to the ISchulz. Mandelbaum fc Padmanabhan! 

l|2010l) results. This fully supports our find ings. 

We further compare our measurements to Ivan Uitert et al.) 
(2011) who measured the GGL signal using the overlapping 
region between the SDSS-R7 and the RCS2 survey. While 
the lens sample is chosen spectroscopically, the shear 
estimates are derived on the deeper RCS2 r'-bands from 
CFHT. Comparing the corresponding luminosity bins we 
observe a good agreement for the excess s urface mass den- 
sity in our measurements and the results of lvan Uitert et al] 
(! 201l j) up to a luminosity of M^/_ab ~ -23 to -24 (cf. their 
fig. 8). 

The dependence of the measured excess surface mass 
density on the environment density is shown in Ap- 
pendix IB2I (see Figs. IB3IB8|I . At this point we then also 
show how well the observations (for either red or blue 
galaxies residing in different environments) are described 
by the best fitting BBS- and NFW halo parameters derived 
for red and blue galaxies later on in this Section. 



4-2.1 L-a-Scaling Based on Fits to AE 

We now fit an isothermal sphere to the AE-profile for the 
previously defined magnitude intervals (see Fig. I26|) . The 
corresponding best fit velocity dispersions are shown in 
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Figure 26. Excess surface mass density AS for individual luminosity bins. The fainter luminosity bins are dominated by blue lens 
galaxies (blue squares, dotted fit-line), outnumbering red galaxies (red triangles, dashed fit-line) in those redshift bins. The red galaxy 
sample shows a significantly higher signal, both due to the higher mass of red galaxies for given luminosity and the denser environment 
red galaxies populate in general. For brighter luminosity bins the fraction of red galaxies increases. 



Fig.[27|as a function of the weighted mean (see eauatioiil62p 
luminosities of galaxies from the respective magnitude 
interval. The black data points show the complete sample 
(all galaxies in all density environments). The results for 
galaxies in high and low density environments are shown 
in magenta and green. To avoid contamination by neigh- 
bouring galaxy haloes, especially in the fainter luminosity 
bins, we only consider AS data points out to ~ 100 



kpc. For bright galaxies the velocity dispersions agree 
irrespective of the galaxy environment density. In high 
density environments faint galaxies show a mildly higher 
signal than their counterparts in low density environments 
(see Fig. [27]), increasing in difference with decreasing 
luminosities. The interpretation for this is that for faint 
galaxies the environment could impact the AE signal 
already on scales smaller than 100 kpc, but more 
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Figure 28. Velocity dispersion ct as a function of luminosity. We determine the velocity dispersion by fitting an SIS out to a scale 
of 100 h^^ kpc to the excess surface mass density AE in separate luminosity bins (see Fig. I26I I. The left panel shows the results for 
the combined lens sample (black circles, solid fit-line) and for red (red triangles, dashed fit-line) and blue (blue squares, dotted fit-line) 
lenses separately. The slopes are (0.24 ± 0.03) and (0.23 ± 0.03) for the red and blue sample. Due to the excess of elliptical galaxies at 
high luminosities and the excess of spiral galaxies at low luminosities this combines to a slope of 0.29 ± 0.02 (black dashed line) for the 
complete lens sample. The right panel shows the corresponding values for a assuming a luminosity evolution according to L oc (1 -I- z). 
The amplitude is increased in comparison to the values without evolution, however, the slope and thus the scaling relation between 
a and luminosity remain almost unaffected (0.30 ± 0.02 for all galaxies, 0.26 ± 0.03 for red and 0.25 ± 0.03 for blue galaxies). The 
parameterisations for the derived scaling relations are shown in Table [Tl 



likely it is the change of fraction of blue and red galaxies: 
in dense environment the red galaxies are relatively more 
abundant and thus increase the AE-amplitude of the 
combined sample. Therefore the scaling relation of the 
velocity dispersion becomes slightly steeper with decreasing 
environment density, being a oc 7;^0;3i±o.o3 density 
environments, a oc /,0;29±o.02 ^^jj density environments 
and finally a oc 7;,o.27±o.02 high density environments. 
With these power law slopes we obtain for a galaxy with 
luminosity L* = 1.6 x 10^" L^/ q, which corresponds to 
an absolute magnitude of M*, ~ —21.7 in the vega system, 
a velocity dispersion of a* = 135 ± 2 km s~^ if all lenses in 
all environments are considered, a* — 141 ± 2 km s~^ for 
all lenses in dense environments and a* = 132 ± 2 km s~^ 
for all lenses in low density environments. 
As next step we analyse the lens sample for red and blue 
lens galaxies separately without further discriminating 
the environment density. Both galaxy types, red and 
blue, show the same scaling behavior, (Trcd oc 7;,0;24±o.03 
and cTbiue oc j;^0-23±o.03 respectively, consiste nt with the 
Faber- Jackson relation l|Faber fc JacksonI Il976l ). as can be 
seen in the left panel of Fig. [2H1 (red and blue dashed lines) . 
However, red galaxies show a significantly higher amplitude 
than their blue counterparts over the whole luminosity 
range. The transition from red SED-dominated galaxies 
at high luminosities to blue SED-dominated galaxies at 
fainter luminosities leads to the steeper scaling relation of 
a oc 7;,0;29±o.o2 £^j. ^j^g complete lens sample (see black solid 



line in the left panel of Fig. [28)) . 

Based on our fit-values for a in the different luminosity 
bins and the obtained scaling relations, we obtain a velocity 
dispersion of a* — 135 ± 2 km s~^ for the combined lens 
sample, a*^^ = 162 ± 2 km s~^ for red lenses and a velocity 
dispersion of a^i^^ = 115 ± 3 km s~^ for blue galaxies with 
luminosity L — L* . 

Until now we have grouped galaxies together that have 
the same rest frame luminosity, ignoring the fact that 
their luminosity evolves with lookback time and thus 
redshift. If we want to group galaxies together that have the 
same present day luminosity we have to account for their 
luminosity evolution from the redshift of observation to 
z = 0. In the following we assume that galaxy luminosities 
evolve as L oc (1 + z). This assum ption is correct for red 
galax ies (see e.g. ISaglia et al. I I2OIOI or also iBernardi et al.l 
|2010D . Blue galaxies evolve more rapidly. For an estimate 
we could take their SED-types and plausible star formation 
histories to obtain their luminosity dimming. This will not 
be extremely precise, and therefore in the following we take 
the evolution of red galaxies as a lower limit. 
Accounting for luminosity evolution increases the veloc- 
ity dispersion for a galaxy with present day reference 
luminosity to a* = 150 ± 2 km s~^ for all galaxies, 
""rod = 1'''3 ± 2 km s~^ for red and fJiuo ~ 123 ± 3 km s~^ 
for blue galaxies. The scaling relations, however, are nearly 
unaffected. We obtain slopes of a oc i0.29±o.02 
galaxies, a oc lO-25±o.o3 ^^^^ ^ ^ ^o.24±o.o3 f^j. galaxies with 
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Figure 27. Velocity dispersion cr as a function of absolute lumi- 
nosity for galaxies in different environments. We determine the 
velocity dispersion by fitting an SIS out to a scale of 100 
kpc to the excess surface mass density AS in separate luminosity 
bins (see Fig. l26|l . We show the results for the combined (red and 
blue galaxies) sample for all environments (black circles, solid fit- 
line), for dense environment (magenta triangles, dashed fit-line) 
and environment with low density (green squares, dotted fit-line). 



red and blue SEDs (see also Table [T|. The closer inspection 
of red and blue data points in both panels of Fig. 1281 
shows that for a luminosity of L ~ 6 — 7 x 10^" Lq 
two red data points are decreased relative to the red 
SED linear fit and that one blue data point is increased 
relative to the blue SED linear fit. This could point to 
a problem in contamination of the red and blue samples 
with blue and red galaxies at this luminosity. We will 
see this feature (that apparently two data points for the 
red galaxies are biased low) in most of the following Figures. 

The velocity dispersion cr^J;^ obtained from the weak 
lensing analysis (out to 100 kpc) describes the circular 
velocity Hdrc.haio = o"haio\/2 of the dar k matter halo 
assum ing an SIS-profile (see equation [7|). iGerhard et al.l 
|20o3) (see their fig. 2) have studied the circular velocity 
curves of local ellipticals with stellar dynamics out to a 
few (^C 3) effectiv e radii. They constrained the anisotropy 
profiles /3(r) (see iBinnev fc Tremaind 1 19871 ) of the stellar 
orbits and obtained that the mean values for {3 are typically 
between 0.2 and 0.4. The detailed dynamical models yield a 
relation between the central stellar velocity dispersion and 
the maximal rotation velocity profile of 



0"sta 



= 0.66 Umax 



(64) 



The radii where these maximal velocities are reached are of 
order of 0.5 times the effective radii. The rotation velocities 
for larger radii (> 7?e) are flat and have values of « 0.9 Um^"- 
If one sets these 'asymptotic values' equal to the halo circu- 



lar velocity we obtain 

l^circ.halo = \/2 a^^a = 0.9 Unix = 0.9 X 1/0.66 (Tsta 



WL 
O" halo 



0.96 a,ta 



(65) 



If one sets the maximal circular velocity equal to the halo 
circular velocity one obtains 



O"halo = i-07 (Tsta 



(66) 



In Fig. [29] we compare how the measured velocity disper- 
sion (Tstar of LRGs comparc with pr edictions from our WL - 
analysis f or red ga l axies, based on lEisenstein" et all (|200lD 
and Gall azzi et al.l (|2006l ). i.e. we add the best-fitting lines 
for the crj^l^-luminosity relation, rescaled with 1/0.96 (in 
magenta) and 1/1.07 (in red). We have added the rela- 
tion between the Us tar and evolution corr ected luminosi- 
ties of SDSS-LRGs faisenstein et al.l I2OOII ) obtained from 
iGallazzT et al.l (|2006l ) as green dashed line. This relation 
is however obtained from fitting a linear relation of veloc- 
ity dispersions vs. absolute magnitude to the overall LRG 
sample. In the lensing analysis we first average the signal 
within some (small) luminosity bin and then study the sig- 
nal. To treat the LRG-galaxies in a similar way we have 
obtained the (Tstai-values from the SDSS data base and es- 
timated the luminosity evolved redshift zero absolute mag- 
nitudes in the r-band (from SED-fits and a luminosity evo- 
lution proportional to 1 -\- z) and obtained the mean stel- 
lar velocity dispersion within equidistant luminosity inter- 
vals. For this we only include galaxies with redshifts be- 
tween 0.05 and 0.3 and with secure velocity dispersion es- 
timates 0.03 < d(Tatar/o"star < 0.1. The results are plot- 
ted with filled black circles, whereas the density contours 
for all considered galaxies are shown in black. We see that 
0.96 (Tstar ^ o-hik) ^ 1-07 (Tstar holds at least for for lumi- 
nosities above 10^" Lq. Therefore the halo velocity is 
between the maximal circular velocity found around 0.5 i?e 
and 90 per cent of this value which equals the velocity of 
galaxies at a few effective radii. This indicates that at least 
for galaxies above this luminosity threshold the halo indeed 
is isothermal on scales out to 100 kpc. 



4-2.2 L-r2oo- and L-K'hm-Scaltng Based on Fits to AE 

We now assume that the co ncentration-mass- r elation 
c oc M-° °®'' (see equation [Ml) of [Duffy et all l|2008l ) holds 
which reduces the NFW-profile to a one-parametric profile. 
We fit NFW-profiles to the AE-profiles (shown in Fig. ^ 
and thus obtain the virial radius for each luminosity inter- 
val. The r2oo vs. luminosity relations (with and without 
luminosity evolution) are shown in Fig. 1301 They imply that 
if the assumed concentration-mass-relation is correct then 
the r2oo vs. luminosity relation cannot be described by a 
single-power law anymore, but instead with double-power 
laws and a break at around L — 10^° Lq. In this 

case, the mean mass-to-light ratio of galaxies within a 
luminosity interval would indeed be minimal at this break 
luminosity luminosity. This is in agreement with results 
from abundance matching (AM) t echniques and som e 



satellite kinematic results (see fig . 1 Dutton et al. [2010l ). 



in particular with the results of [More et al.l (20l3) (see 



their fig. 5) who also obtained a change of slope for the red 
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Figure 30. Virial radius r200 as a function of absolute luminosity. The left panel shows the result without, the right panel with luminosity 
evolution L oc (1 + z). The red triangles and dashed fit-lines denote red galaxies, the blue squares and dotted fit-lines blue galaxies and 
the black cirlces and solid fit-lines all galaxies. We only use galaxies with L > 10'^'^ Lq for the determination of the scaling 

relation. For the combined lens sample the r2oo scales with i0.39±0.04 ignoring, and with £,0-37±0.04 including luminosity evolution. The 
parameterisations for the derived scaling relations are shown in Table [Tl 
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Figure 31. Circular velocity V200 as a function of absolute luminosity for red (left panel) and blue galaxies (right panel). Analogously to 
the fit of the r2oo we onl y use data points w ith L > 10^'' Lp, for the determination of the scaling relation fit. In the right panel we 
have added the results of lReves et alj l|201ll ) and the model of lDutton et al. (201 o|), after tr ansl ating their stellar mas ses to luminosities 
fo r blue galaxies (see t ext). In the left panel we have added the results of lOutton et al.l | |201C|) and lGallazzT et al. 1 1I2OO6I I using the relation 
of lDutton et al.l l|201Ch (see text). 
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Figure 32. Virial mass M200 as a function of absolute luminosity. The left panel shows the result without, the right panel with luminosity 
evolution L oc (1 + z). The red triangles and dashed fit-lines denote red galaxies, the blue squares and dotted fit-lines blue galaxies and 
the black circles and solid fit-lines all galaxies. We see as expected the same scaling behavior as for the T200 (see Fig. I30| l. Analogously 
to the fits of the r2oo we only use data points with L > 10^" Lq for the determination of the scaling relation. For the complete 
lens sample the M200 s cales with £,121±0.10 ignoring, and with /,i i2±0.ii including luminosity evolution. We included th e resul ts from 
Ivan Uitert et al. I l l201ll'l for a red SED-dominated lens sample in the right panel (green crosses) and from iHoekstra et al.l ||2005| ) for an 
isolated and thus blue SED-dominated galaxy sample (open magenta circles). The parameterisations for the derived scaling relations are 
shown in Table [T] 



Table 1. 



Without luminosity evolution, L*, = 1.6 X 10^" h ^ ^r',© 



typo 


(T* [km s 




r|(,(, [h 1 kpc] 


'?''200 


MSoo [10" h-^ 


Mq] 


VM200 


All 


135 ± 2 


0.29 ±0.02 


146 ± 2 


0.39 ± 0.03 


11.1 ± 0.4 




1.21 ±0.10 


Red 


162 ± 2 


0.24 ±0.03 


177 ±3 


0.33 ± 0.04 


18.6 ± 0.8 




1.05 ±0.12 


Blue 


115 ± 3 


0.23 ±0.03 


120 ± 2 


0.36 ± 0.07 


5.8 ±0.5 




1.14 ±0.20 






With luminosity evolution, L*, = 


1.6 X 10^" h- 


ir',0 X (1 + z) 






type 


(T* [km s"-' 


»7<T 


^200 ['^ ^ kpc] 




M|oo [10" h'^ 


Mq] 


VM200 


All 


150 ±2 


0.29 ±0.02 


170 ±2 


0.37 ±0.04 


17.0 ±0.6 




1.12 ±0.11 


Red 


173 ±2 


0.25 ±0.03 


198 ±3 


0.38 ± 0.04 


26.1 ± 1.1 




1.17±0.13 


Blue 


123 ±3 


0.24 ±0.03 


133 ±3 


0.40 ± 0.08 


8.7 ±0.6 




1.37 ±0.25 



Exponents for the scaling relations of the velocity dispersion tr, assuming an SIS and for the r200 and Af200j assuming an NFW profile 
without and with luminosity evolution. The SIS fits have been extracted from all all luminosity bins, the fits for the NFW profiles only 

include luminosities brighter than L = 10^" Lq. 



galaxies' M200 vs. luminosity relation at a luminosity of 
about 10^" Lq. We would like to point out however, 
that the result in Fig. |30] only holds if the concentration is 
only weakly changing with virial mass. One could instead 
approximately reconcile a single-power law r2oo-luminosity 
relation if one required the concentration to rise steeply for 
luminosities smaller than 10^^° Lq. We will investigate 
these two alternatives more in Section 14.3.31 Because of 
the apparently broken r2oo-luminosity scaling relation we 
measure the power law slope only for galaxies brighter than 



10^° h-^ Lq. We obtain power laws of r^^o oc lO-33±o.04 
and r^ii^" oc L°-^^^° °'^ for red and blue galaxies with- 
out luminosity evolution and of r2oo oc 7;,0 38±o.04 
r^oo" oc L° '"'±°"** for luminosities evolving with (1 + z). 
If galaxies are not separated into blue and red SED types 
we obtain (for the combined sample) r2oo oc 7;,0-39±0 03 
and r-200 oc /^0-37±o.04 ;gjjQ]-;jjg luminosity evolution and 
assuming a (1 + 2) scaling. As before the steeper scaling 
is due to the fact that the amplitudes for the r2oo vs. 
luminosity scalings are different for red and blue galaxies 
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Figure 29. Velocity dispersion for our red galaxy sample, 
rescaled by 1/0.96 (in magenta) and 1/1.07 (in red) as a func- 
ti on of absolute l umino sity. We compare our values to the ones 
of lGallazzi et al.l ||2006|) (ereen da s hed fi t-line) and a SDSS-LRG 
sample based on lEisenstein et al ] l|200lh (black circles and con- 
tours), only considering LRGs with redshifts of 0.05 < z < 0.3 
and velocity dispersion uncertainties of 0.03 < do-star/cstar < 0.1. 



and the fractional mix of red and blue galaxies changes as 
a function of luminosity. 

In Fig. [31] we translated the result for r2oo to the virial ve- 
locity 11200 using equation (|26p . In the right panel of Fig. [31] 
we show V200 vs. luminosity for our blue sample (blue data 
points) and t he power law fit fo r L > 10^'' Lq (blue 
dotted line). [Reyes et al.l (|201ll ) have measured D200 for 
SDSS disc galaxies as function of stellar mass. We try to 
compare their result to ours by translating their stellar mass 
estimate (ba ck) to luminosity. It seems that for local disc 
galaxies (the [Reves et al.|[201ll disc galaxies have redshifts 
between 0.02 and 0.1) an on average mass-to-light ratio of 
Mstai/Lr — 1 Mq/Lq is a go od description. This can be 
seen on one hand in fig. 1 of Ivan Uitert et al.l l|201ll ) by 
comparing their blue histograms on the vertical to the hor- 
izontal axis showing the luminosity distribution and stellar 
mass distributio n of blue SDSS-g alaxies. This estimate is in 
agreement with [Bell et al.l l|2003l ). if one takes into account 
that our local (see Fig. I21|l galaxies have a. {g — r)-colour 
of approximately 0.3 — 0.4 at the b right end (which are 



the galaxies in common with iReves et al. 201Jj}. The same 
result is obtained from iKauffmann et al.l (j 20031 ) . fig. 14, 
upper right panel, if one takes into account that our local 
blue galaxies have mostly absolute magnitudes fainter than 
M^ / = —21. Fo r the three luminosity intervals provided 
by [Reves et al.l (|201ll ) their data points (translated to 
luminosity) agree well with ours (see Fig. 1311 right panel). 
We have a larger dynamical range and can extend our 
analysis down to to a few ti mes 10^ Lp,. In an an alogous 
way we have translated the [Putton et al.l l|2010l ) model 



for the 11200-stellar mass relation to the W2oo-luminosity 
relation. It agrees very well with our result, but might have 
a slightly shallower slope. 

For the red galaxies we have translated our luminosi- 
ties into stellar mass estimates, because in this case 
we are confident that we can do so for the abso- 
lute magnitude and redshift range considered. We use 
loKin(Mitar) = 1.09 3 logj^Q Lr — 0.57 3 (which was used b y 
[Putton et all I2OI0I and derived from [Gallazzi et al.lliooel ). 
and insert luminosity evolution corrected luminosities. Our 
results for 1120 are shown i n red in Fig. 1311 together with 
the model of [Putton et al.l (|2010l ) . They are the same to 
a remarkable level. Only the results for the second and 
third brightest luminosity interval lie below, for a reason we 
speculated about already. O n top we have a dded the result 
for Vapt as obtained from [Gallazzi et al.l (|2006[ ) velocity 
dis persion vs. l u minos ity relation, using the prefactors 
of [Putton et al.l l|2010l ) for the relation between velocity 
dispersion and rotation velocit y. (Using our translation 
factor the [Gallazzi et al.l [20061 curve would be slightly 
below). We concluded that for luminosities between 10^" 
and 6 x 10^° Lq the mass density profile of ellipticals is 
not only isothermal out to 100 kpc (as shown before), 
but out to the virial radius. For higher luminosities the 
virial velocity exceeds the optical velocity. 

Finally we translate our r2oo vs. luminosity results 
into virial masses and show results with and without 
luminosity evolution correction in the left and right panels 
of Fig. 1321 For the power law fits (added as red dashed 
and blue dotted lines) we have again used only galaxies 
with L > 10^" Lq. For the combined sample we obtain 
M200 oc _Li-2i±o io and M200 oc L^-'^^±o.ii 
out and with luminosity ev olution correctio n . This scaling 
agrees with the resuhs of iGuzik fc SeliakI (|2002[ ) within 
their larger un certainties (M o c L^^^^"'^''). We have added 
the results of iHoekstra et"aLl (|2005l ) as magenta points, 
which agree well with our blue sample. This ag reement 
appears reasonable since the iHoekstra et all (120051 ) sample 
contains isolated galaxies which mostly consist of blue 
galaxies. In addition we have considered t he ex cess surface 
mass density profile of Ivan Uitert et al.] ([201 ll ) (see their 
fig. 8), and translated them into virial mass estimates in 
the same way as we did for our work. These estimates are 
shown as green points. They agree well wit h our red sample 
result s, which again is reasonable since the Ivan Uitert et al.l 
(|201ll ) sample is dominated by red galaxies. All results 
obtained for r2oo and M200 are summarised in Table 2. 
At last we translate our A/20o vs. luminosity rela- 
tion from the right panel of Fig. [32] into the M200 
vs. stellar mass relation (MSR) using the relation 
log^Q(Mstar) = 1.093 log^Q Lr — 0.573 as above. The result 
is shown in Fig. 1331 The virial-to-stellar mass ratio (shown 
as red points) is almost constant (at 100) for a decade in 
stellar mass (10^" to 10^^ h-^ Mq), and increases for lower 
stellar masses. In the mass range of 10^° to lO'^'^ Ma\ 



our result precisely agrees with the [Putton et all ([2010D 
model shown as black solid curve. At the high stellar mass 
end the MSR seems to increase (if at all) only slowly 
with stellar mass. This satur a tion is in agreement with 
the results of Ivan Uitert et all (|201ll ) ( green points, taken 
from their fig. 14, and converting their stellar masses to 
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Figure 33. Stellar Mass vs. A'/20o/A'^star ratio for red galaxies 
converted to 2 = 0. The red triangles denote our r ed galax - 
ies. We have add ed the results of Mandelbaum ct al.l l|2006d ) 
(open tri angles'). iMandelbaum. Seliak fc Hirata t200a) (open 
squares), I Schulz, "MandeI baun7fc|Padma^^ ||2010| ) (black 

crosses) and DuttoiT^ alj 1I2OI0I ) (black solid line), see also 
fig. 1 in iDutton et al. I l l201Cf) . We further include the results of 
Ivan Uitert et al. I | |2011D from their fig. 14 as green open circles. 



Ho = 100 km s"^ Mpc"\ as in this Figure the stellar 
masses are given for Hq — 70 km s~^ Mpc~^ and the virial 
masses are given for Ho — 100 km s~^ Mpc~^ according to 
van Uitert, private communication) which also saturates at 
a value of about 100 to 150. The Ivan Uitert et all (|201ll ) 
point s for low stellar ma s ses are however even below the 
early IMandelbaum et al.l ll2006d) res ults and seem very 
low. Since the Ivan Uitert et al.l ( 20111 ) M200 vs. luminosity 
relation derived by us from their AE results agree well with 
ours, the dilference can only be due to a different rela tion 
for the stellar masses (Note, that Ivan Uitert et al.1l201ll aim 
to add up the total stellar mass, i.e. not only that of the 
central galaxy but also that of its satellites). 



4.3 Maximum likelihood analysis 

In the previous subsection we have analysed the lens- 
ing signal on small scales, ignoring the impact of further 
galaxy haloes. We now analyse the signal out to scales 
R ^2 Mpc in the case of BBS. Since the integrated 
mass for the NFW profile diverges for infinite radii and the 
integrated mass value within a radius of 1 Mpc already 
exceeds the total BBS mass, assuming reference halo pa- 
rameters of a* ~ 130 km s~^, s* ~ 200 kpc, c* ~ 6 
and r2oo ~ 130 kpc, we limit the NFW maximum like- 
lihood analysis to a maximum distance of R ^ 400 h^^ kpc. 
A projected separation oi R = 400 kpc corresponds to 
~ 1 r2oo for bright and ~ 5 r2oo for faint galaxies. Since at 
these scales the original NFW profile hardly differs from a 



truncated NFW profile (see lBaltz. Marshall fc Ogurill2009l ) 
the not-finite total NFW masses do not affect our results. 
At scales of 400 kpc the signal of other haloes becomes 
comparable to the haloes of interest. This means, we have to 
model all the haloes simultaneously. We do this in a maxi- 
mum likelihood analysis (see Section r2.5p . where we measure 
the (two) parameters of a reference galaxy (with a reference 
luminosity) and describe haloes of galaxies with different lu- 
minosities by scaling relations. As reference luminosity we 
choose Lr' = 1.6 x 10^ 



h ^ Li,'_0 for the rest of this work. 



4.3.1 Truncated isothermal sphere (BBS) 



In Section 14.2.1! we had analysed the combined sample 
fitting an SIS to the excess surface mass density AS 
and obtained a scaling relation of ct oc j,0-29±o.02 
combined sample in Section 14.2.11 As a cross-check we 
now analyse the luminosity vs. velocity dispersion scaling 
with a maximum likelihood analysis. We use an SIS and 
analyse foreground-background pairs only to a distance of 
200 kpc, since in this region the difference between the 
singular and the truncated isothermal sphere is negligible. 
In this way all (visible) nearby haloes are accounted for, 
in particular the environment of satellite galaxies sitting in 
more massive haloes traced by another galaxy. This analysis 
yields an exponent of ijo- = 0.31 ±0.02 (see black contours in 
Fig. I34|) . in agreement with the results from Section [4.2.11 
The inclusion of luminosity evolution with L oc (1 -I- z) 
increases the values for ct* (see also Section I4.2.ip . but 
basically does not change the slope of the scaling relation, 
as the best fit values (triangles in Fig. I34p show. Therefore 
we now fix this scaling to ct oc L"'^ for the combined 
sample. In Section [4.2.21 we have found AI200 oc L^'^ for the 
combined sample. Although NFW and BBS masses are not 
the same, we expect similar lu minosity scalings to hold. 
We th erefore, in agreement with iHoekstra. Yee Sz GladdersI 
(|2004l ). assume Mbbs oc L^'^ to hold. This then (using 
equation (31} implies 77s = 0.6. 

With these assumptions and if we ignore luminosity 
evolution we obtain for the reference L*-galaxy a veloc- 
ity dispersion of ct* = 131^2 km and a truncation 
radius of s* — 184:'^\'[ h~^ kpc, implying a BBS mass of 
-^totai.BBS = 2.4l°;2 X 10^^ h~^ M© (measuring the haloes 
out to scales of 2 h~^ Mpc). The results for the case of 
luminosity evolution are given in Tabled 
The left panel of Fig. [21] showed that the luminosity vs. ct 
scaling of the combined lens samples differs from that of 
pure red or blue galaxy samples, as it additionally traces 
the transition from massive red SED-dominated galaxies 
to low-mass blue SED-dominated galaxies. We therefore 
separately investigate the scaling relations for red and blue 
galaxies. In the analysis of red and blue galaxies, only 
for one galaxy type at a time the best-fitting parameters 
are examined, while for the other galaxy type the lensing 
signal is calculated with adjusted values for the fiducial 
quantities, as e.g. the velocity dispersion ct and the scaling 
parameter r]cr. I.e. we assume fixed values ct* and r]a for 
the blue galaxies and run a maximum likelihood for the 
corresponding values for red galaxies, and vice verse. These 
results then are each fed into a further iteration for the 
maximum likelihood analysis, repeating the procedure 
until the calculations converge and do not deliver results 
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Figure 34. Constraints on velocity dispersion a* and the velocity dispersion vs. luminosity scaling index 77^ for a fiducial galaxy with 
luminosity L*, = 1.6 X 10^^ L^i q. Blue galaxies are shown in blue in the left, red galaxies in red in the right and the combined 
galaxy sample in black in the middle. The contours indicate the 68.3 per cent (solid), the 95.4 per cent (dashed) and the 99.7 per cent 
(dotted) confidence levels for both parameters combined, the crosses indicate the best-fit values, the filled triangles indicate the best-fit 
values assuming luminosity evolution with (1 + z). 



deviating from the previous iteration step. We obtain values 
of rj'J"^ = 0.27t^;J!^ for red and 77^'"'= = 0.27to:ol for blue 
lenses (see Fig. I34|l . 

We now fix the slope of the red and blue samples to 
V^'^ ~ '7ct'"'^ ~ 0.25 and measure the velocity dispersion and 
truncation radius of the red and blue reference galaxies. 
For this we keep the assumption regarding size vs. velocity 
dispersion scaling of s oc ct^. Using -qs = 0.5 in equation (|3ip 
then gives a M/L-ratio independent from luminosity for 
both the red and blue galaxy samples. 

The analysis of the red lens sample yields a velocity 
dispersion of a*^^ = I49I3 '^^'^ ^ truncation 

radius s*^^ = 887^37 kpc, implying a total mass of 
^tot^BBS = 5-5lg;g X 10^^ h-^ Mq. For the blue lens sam- 
ple we obtain a velocity dispersion of a^i^^ — HSlg km s^^ 
and a truncation radius s^iuo ~ 841^4 kpc and a total 
mass of Mj*^^''^'j3Bg = S.elai x 10^^ Mq These results 
are also visualised in the upper left panel of Fig. 1851 
We measure the same quantities for galaxies in high and 
low density environments. It can be seen in the lower 
panels of Fig. [35] that the environment hardly influences 
the velocity dispersion of the galaxies. This means that 



the central matter density of galajcies mostly depends on 
their luminosities and hardly on their environment. The 
truncation radius however, increases significantly with the 
density of the environment (see Fig. as well) for galaxies 
of both SED-types. The best-fitting values for the different 
samples are shown in Tabled 

In general dense environments make it more difficult to 
measure the truncation radius s, because at distances where 
the 'end' of the halo becomes measurable many other haloes 
then contribute to the combined signal. Further we would 
like to point out that the interpretation of the BBS-mass 
(and thus the size s) is that it gives the total mass for a 
halo hosting a galaxy with a central density profile with 
amplitude a. Since in dense environments there is also dark 
matter on scales of groups and clusters the mean mass (and 
thus the truncation radius) associated to galaxies grows. 

4-3.2 Universal density profile (NFW) 

In Section 14.2.21 we carried out an NFW fit to the AE- 
profiles on scales R ^ 100 kpc and obtained a scal- 
ing compatible with r2oo oc L"'* for the combined galaxy 
sample and r2oo oc L^^^ for the red and blue galaxy sam- 
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Figure 35. Maximum likelihood results for the BBS profile. The truncation radius s is shown on the x-axis, the velocity dispersion a* 
is shown on the y-axis, the contours show the 68.3 per cent, the 95.4 per cent and the 99.7 per cent confidence levels. In the upper row 
the results for different galaxy types are compared for one environment in each panel. Black contours show the combined galaxy lens 
samples, red contours the red and blue contours show the blue galaxy lens samples. The left panel shows the galaxies in all environments, 
the middle one shows the galaxies in high density environment and the right panel shows the galaxies in low density environment. In the 
lower row we compare the results for identical galaxy types in different environments, black is all environments, magenta is high density 
environment and green is low density environment. The left panel shows all galaxy types, the middle one shows red galaxies and the right 
panel shows blue galaxies. We see that red galaxies in general have larger velocity dispersions and halo sizes than blue galaxies, leading 
to higher masses. This holds for all environments. Looking at the dependence of the individual galaxy profiles with respect to their 
environmental conditions no galaxy type shows significant evolution in velocity dispersion, but the truncation radius clearly increases 
with environment density. 



pies. We now fix these slopes, i.e. we assume the mass- 
luminosity relation to scale like this for the overall lumi- 
nosity range (single power law) a nd we also a ssum e the 
concentration-luminosity relation of iDuffy et al.l (|2008| ) (see 
equation I36p to hold. The maximum likelihood analysis for 
the combined sample yields c* — 6.41q ^ and a virial radius 
of fJoo = 133l2 kpc. This translates into a virial mass 
of M200 = 7.6l°:^ X 10^1 ft-i Mq. 

As next step we measure the parameters for the SED 
types separately. The procedure is the same as for the 
BBS profile: we fix one galaxy type with adjusted fidu- 
cial parameter values (assuming M/L=constant) and ob- 
tain the best-fitting parameters for the other galaxy 
type, iterating the analyses until convergence. For red 



SED lenses we measure c*''^'"^ = 6.4lQg and a virial 
radius of rjoo"'* — I6OI4 kpc, leading to a virial 
mass of M*oo°'* = l.2to\ x 10^^ /i"^ Mq. For blue galax- 
ies, we measure c*''^'"'^ — 7.0ljg and a virial radius 
'"260'"° ~ f fSlg kpc, resulting in a virial mass of 
^^260'"'' = 5.0lo;6 X 10" h"^ Mq. The virial radu of the red 
galaxies are significantly larger than for the blue galaxies, 
leading to the significant mass difference, as already mea- 
sured for the BBS profile. This also can be seen in the up- 
per row of Fig. 1361 For the (same) fiducial luminosity the 
halo concentration thus is similar for red and blue galax- 
ies. But since the virial radii and thus the virial masses for 
the red galaxies are significantly larger than for blue galax- 
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Figure 36. Maximum likelihood results for the NFW profile. The concentration parameter c* is shown on the x-axis, the rjQQ is shown 
on the y-axis. The contours show the 68.3 per cent, the 95.4 per cent and the 99.7 per cent confidence levels. In the upper row different 
galaxy types are compared for the one environment in each panel: black contours show the combined galaxy lens sample, red contours 
the red and blue contours show the blue galaxy lens sample. The left panel shows the galaxies in all environments, the middle one shows 
the galaxies in high density environment and the right panel shows the galaxies in low density environment. In the lower row we compare 
the results for identical galaxy types in different environments, black is all environments, magenta is high density environment and green 
is low density environment. The left panel shows all galaxy types, the middle one shows red galaxies and the right panel shows blue 
galaxies. The concentration shows no significant difference comparing the different galaxy types at the reference luminosity, but the r2oo 
is significantly higher for red galaxies than for blue galaxies, also increasing with environmental density for all SED types. 



ies (of the same luminosity), the concentration is higher for 
massive red galaxies, if galaxies with the same masses are 
compared. If we scale the mass of a blue galaxy to a value 
of MjYo"'^ = 12.4 X 10" h-^ Mq, we obtain a concentration 
parameter of c'''"" = S.IIq'^ for this galaxy. For less massive 

galaxies (Af200 < Afjoo'"'') this relations inverts. The results 
of the maximum likelihood analyses are summarised in Ta- 
ble H 

Investigating the dependence of NFW haloes on the envi- 
ronment (see also the lower panels of Fig. I36p shows that 
galaxies in high density environments hardly differ in con- 
centration and only have a slightly higher virial radius than 
galaxies in average density environments, but that galaxies 
in low density environments have a significantly lower virial 
radius and also show a tendency of having a higher con- 
centration (especially blue lens galaxies) than in denser en- 



vironments. The best-fitting values for the NFW-likelihood 
analyses are shown in Table [5] 

We finally check the consistency of our BBS and NFW 
analyses by comparing the corresponding (approximate) 
virial masses. We calculate the halo masses based on 
the profile parameters within the r2oo (obtained from 
the NFW fitting). For the maximum lens sample the 
BBS profile yields MBBs(r20o) = 9.2lo9 x lO^^fe"^ Mq, 
which is only slightly higher than the corresponding 
Af2oo from the NFW analysis. For the total lens sam- 
ple in high and low density environment we mea- 
sure a mass of MBBs(r|oo) = 10.lli;4 x 10"/i~^ Mq and 
AfBBs(''^20o) = S.SlJ^'g X lO^/i""*^ Mq, respectively, which is 
still slightly higher than the corresponding NFW values, but 
agrees well within the uncertainties. 

The same comparison for the red and blue galaxy lens sam- 
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Table 2. 



Without luminosity evolution, L*, = 1.6 X 10""' h ^ L^/ q 



type 


density 


[km s^^] 


s* 

[/i-i kpc] 


c* 


^"200 
kpc] 


^^total BBS 

[10" h-'i Mq] 


[10" /l-l Mq] 


[10" /l-l Mq] 


A^totalBBs/-^* 

[/I M0/L,,,q] 


all 

all 
all 


all 
high 
low 


131+3 

13lt^ 


184111 


6.41°;? 

6.4ti;[! 

9.4t2.4 


133^2 
137+1 


23.2t2'5 
09 9+4.5 

19 1+2.6 


7-6i;jj 

q+0.5 

c; q+0.5 


„ ^-1-1 1 

9-2lS;^ 

10.1111 
6.8111 


178lig 

248+3? 
93l^« 


red 
red 
red 


all 
high 
low 


149^3 
150t^ 
144t| 


337l« 
245tg 


fi 4+0-'' 
"■^-0.8 

fi 0+^° 


160+4 
1671^ 
146^5 


54.9tr.^2 

07 q+12.3 

"^-10.5 


14.111:^ 

9.4ll;0 


15.5+?-? 
16.81^4 
12.81^-° 


422l^J 

589111^8 
2871^^ 


blue 
blue 
blue 


all 
high 
low 


1141^ 
126l| 


1071^2 
401- 


7.011-1 
7.0i3- 


115+^ 
105l« 


8.6+ii 

io.2i:° 
4.7ti:? 


0+0.5 
^•"-0.7 

q 0+0.6 


s.itll 
5.311;? 

3.6111 


6611? 
78lg 
3611^ 








With luminosity evolution, L*, 


= 1.6 X 10"' /l-2 


L^/^Q X (1 + z) 






type 


density 


cr* 
[km s^^] 


s* 

[h-i kpc] 


c* 


'"200 
[/i-i kpc] 


-'^'^total.BBS 

[10" /i-i M0] 


-'^-^200 
[10" /l-l Mq] 


.'^-^BBs('"20o) 
[10" h-^ Mq] 


.'^■^total.BBs/-^* 
[h Mq/L,,,q] 


all 
red 
blue 


all 
all 
all 


144^^ 
161^3 
1261^ 


2531^3 
108ii? 


fi 2+°-* 
"•^-0.7 

fi 2+1-^ 
"•^-1.4 


158+2 
183t* 
1351^ 


38.5j;^;l 

12.611:2 


12 7+°-^ 
18.6+1-2 

8.oi?:j 


13.711-1 

20.9121 
7 9+1.8 
'•^-1.5 


296+32 
6051^3 

971^3 



Best-fitting parameter values from the maximum likelihood analyses for the different galaxy SEDs and environmental densities, 
velocity dispersion a" and truncation radius s* for a BBS profile, concentration parameter c* and virial radius rJgQ for an NFW profile, 
the corresponding total mass (BBS) and virial mass (NFW) and the total M/L-ratio based on the total BBS mass. The upper table 
shows the values assuming no luminosity evolution with redshift, the lower table shows the values assuming luminosity evolution with 

L oc (1 + z). 



pies gives analogous results, as can be seen in Table (2] The 
same is true for results obtained assuming L oc {1 + z), 
see Table [2] In general the BBS profiles give slightly higher 
values for the masses t han the NFW profi l es. Th is is con- 
sistent with the results IWright fc BrainerdI (|200Cll ) who ob- 
tained that fitting isothermal profiles to NFW profiles leads 
to higher mass estimates for the same given gravitational 
shear within the virial radius. 

In Fig. [37] we show how well the AE(i?) predictions based 
on the best fit model parameters for the I/*-galaxies and 
the assumed scaling relations fit to the observed AE-profile 
(black). Here the red symbols show the AE-profiles based on 
the BBS simulation whereas the green symbols show the case 
for the NFW simulation. The dashed red and green curves 
show the analytic BBS and NFW single haloes scaled to 
correspond to the considered luminosity bins. 



4-3.3 Scaling relations for the M/L ratio and the 
concentration parameter c 

Until now in our analyses we assumed a scaling relation of 
Mbbs/L oc L"-^ (see eauationl28[) . Given the stability of the 
L — (T-relation we now keep the slope ria fixed and addition- 
ally fix the velocity dispersion to a value of cr* — 131 km s~^ 
for an I/*-galaxy and run a BBS maximum likelihood anal- 
ysis to measure the scaling index r]s for the truncation ra- 
dius s. We find ri, = 0.521°-?^ for s (see black contours 
in Fig. I38|) . As shown in equation (|16|) the total BBS-mass 
scales as M oc a^s, from one obtains the scaling relation for 



the mass-to-light ratio: 

Vm/l = 2ri„ + ris - 1 , (67) 

(see also equation (31} . Combined with the scaling behavior 
of the velocity dispersion rjcr = 0.30±0.02 we obtain a scaling 
of the mass-to-light ratio of 

(M/L) oc l''-i2_o 11 ^gg^ 
for the combined sample. 

The calculation of the scaling relation of s for individual 
galaxy types is challenging and requires a large galaxy 
sample to be statistically significant. This holds especially 
for galaxies in dense environments such as red galaxies, 
because in this case the signal of the outer halo in the 
AE(i?)-profile is superposed by the imprint of neighbouring 
haloes. Since the truncation radius itself is already difficult 
to measure, it is hardly possible to constrain the scaling 
relation unless one considers a very large sample of galaxies 
which are approximately undisturbed by neighbouring 
haloes. We nevertheless try to analyse the scaling relation 
for s for red and blue galaxies separately and show the 
results in Fig. 1381 We see that the scaling relation for the 
blue galaxy samples follows our expectations based on the 
result that »7s'"° = 0.42lg-lJ, which corresponds to a scaling 

relation of M/L oc /."''f/J^ with ri^/)l = -0.08l!l:li (which 
is within the errors consistent to a constant mass-to-light 
ratio). However, considering the red galaxy sample, the 
result of 77^°'^ = — O.I2IQ-I3 is surprising and difficult to 
understand at the first glance. To determine the reason 
we now further investigate our red and blue lens samples 
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Figure 37. Excess surface mass density AS for different luminosity bins. Black circles are the observational data points, red triangles 
result from a BBS simulation, green squares stand for the NFW simulation, based on maximum likelihood analyses (see Section 12.51 1. 
The exact procedure for the simulation is described in Appendix |B] The dashed red and dotted green lines show the analytic single dark 
matter haloes scaled to corresponding luminosity for the BBS and NFW profile, respectively. 



as a function of their rest frame luminosity, splitting the 
samples into bright [M^i < —21) and faint lens samples 
{Mr-i ^ -21). As we see in Fig. [39] the resuHs for the 
blue galaxies remain unchanged, both luminosity samples 
showing the same scaling behavior, i.e. »?s^'bright = 0-50lo;42 
for bright compared to ri^^^Snt ~ 0.45to'io for faint blue 
galaxies, roughly corresponding to luminosity-independent 



mass-to-light-ratio. For the red sample we find that the 
deviating scaling behavior seen in Fig. [3H1 has its origin in 
the faint part of our red lens sample. While the bright red 
galaxies in principle scale as expected (T?s|bright = O.SSlg 21' 
corresponding to »7M/i^bright — ~0-^'^-o^2j consistent with 
a constant mass-to- light ratio), the scaling behavior for the 
fainter fraction is inverted (r^s^ftint ~ '"^■■^^ -o.it)- Since the 
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V, = 1.6x1010 h-2 L . 

r r. 




Figure 38. Scaling relation for the BBS truncation radius s with luminosity. The blue contours on the left show the constraints for blue 
galaxies, the red contours on the right show the constraints for red galaxies and the black contours in the middle show the combined 
galaxy sample. The contours indicate the 68.3 per cent (solid), the 95.4 per cent (dashed) and the 99.7 per cent (dotted) confidence 
levels, the cross indicates the best-fit value. 



velocity dispersion is decreasing with decreasing luminosity 
this implies that the halo masses are only slowly decreasing 
with luminosity. This is in ag reement w it h the halo mass 
vs. stellar mass relation from IGuo et all (|2010l ) (based on 
abundance matching) . Their fig. 2 (upper panel) shows that 
the halo mass decreases only slowly with stellar mass once 
the stellar mass is below ~ 10^'' Mq. To exclude that 
environment structure causes this rise of the truncation 
radius with decreasing luminosity we have repeated this 
analysis for the low density environment galaxy sample. We 
have summarised our results in Fig. 1401 the M/L- ratio for 
red galaxies steeply decreases for increasing luminosities for 
galaxies fainter than L ~ 10^° L® and then turns into a 
further, but shallow decrease for more massive galaxies. In 
contrast, for blue galaxies the M/L-ratio is approximately 
constant for the whole luminosity range considered. The 
confidence intervals shown are based on the uncertainties 
of the halo parameters (for the best fit luminosity scaling 
parameters) and do not include the uncertainties for 
the luminosity scaling parameters themselves. Within its 
uncertainties the slope of the scaling relation is consis- 



tent with values slightly larger than rjs = 0.5, which would 
correspond to a slight increase of M/L-ratio with luminosity. 

Further we consider the scaling relation for the con- 
centration parameter c of the NFW profile with luminosity. 
We perform a maximum likelihood analysis with fixed 
fiducial virial radius rloo, obtained from the previous 
NFW analysis, using the concentration parameter c and its 
scaling exponent rjc as fit parameters. For the combined 

-0 07+"--'-' 

lens sample we find a scaling relation ofcocL ' -on 
(see the black contours in Fig. I4ip . For a scaling of the 

-0.06l-°™ 
-"200 

„ „^ + 0.10 

r0.12 

agreement with the results of Duffy et al.l ('2008'). We also 
investigate possible differences in the scaling behavior for 
red and blue lenses, running separate examinations, and 
repeat the calculations with the results from the previous 
iteration until convergence, as described in the previous 
subsections. Indeed the scaling behavior differs. While the 
concentration parameter of red galaxies hardly depends on 



M2oo/i-ratio with this leads to c oc M^^^ -0.09^ 

06+°'^° 

M200/L oc L°-^^ to c oc M,n n both values being in 



Dark matter halo properties from GGL 43 



_ S 







-1 - 



2 - 



1 - 



- 



"1 — I — I — r 

77„ = 



"1 — I — I — I — r 

= 0.50+0- 



L;. = 1.6x1010 h-2 L^.Q 



"1 — r 



"1 — I — I — r 



"1 — I — r 



"1 — I — I — I — □ 



41 
0.42 



79+22 
' ^ -17 




h-i kpc s* = 291!42 h-i kpd 




- s = 



4f^+o.i2 

^•^^-0.10 

83!;9 h-i kpc 



s = 



0.38!o.{4 
148!^8 h-' kpd 




J L 



J I I L 







100 



200 300 
s* [h~i kpc] 



400 



500 



Figure 39. Scaling relations for the BBS truncation radius s distinguishing between bright and faint foreground lenses. Blue galaxies 
are shown as blue contours in the left side of both panels, red galaxies in the right in red contours. The upper panel shows the scaling 
behavior for bright galaxies with M^i < —21, the lower panel the faint lens sample with M^/ —21. We see that the results for bright 
(r;^ = 0.50l°;42) and faint = 0.45l°-:to) blue galaxies are consistent with each other, following the predicted scaling behavior of 
s oc L"-^. We also see that for bright red galaxies the scaling behavior (r]s = O.SS^q'j^) now also agrees with our expectations within the 
uncertainties, while the deviating scaling behavior of the complete red lens sample is caused by the faint lens fraction, as we see in the 
lower panel {rjs = -0.38^0:13). 



the luminosity {c(x L 



'O Mto.la), the best fit for the slope of 
the c-L-relation of blue galaxies is negative (c cx L -0.26 ^ 
see Fig. I4ip . However, the uncertainties for 7?^"'"° are quite 
large and thus the confidence intervals for ril°'^ and ?7c'"° 
overlap. 



The concentration-luminosity relation can be easily trans- 
lated into a concentration-mass relation, for a given 
luminosity-mass relation. We assume M200 oc L^'^ , but in 
any case the exact value for the scaling exponent is not 
important, as the relation only weakly depends on it and 
hardly changes if e.g. M200 oc i is assumed. Using the 
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Figure 40. Truncation radius s (left panel) and M/L-ratio (right panel), calculated for red (in red) and blue galaxies (in blue) according 
the scaling relation from Fig. 1391 The solid lines denote the best fit, the dashed lines denote the 68.3 per cent-confidence levels (based on 
the uncertainties of the halo parameters and not the scaling relation). The negative s-slope for faint and positive s-slope for blue galaxies 
lead to a strong decrease in the M/L-ratio for faint red galaxies and a weak further decrease for bright red galaxies, while the universally 
positive s-slope for blue galaxies lead to an approximately luminosity-independent M/L-ratio. Within its uncertainties the slope for the 
brighter galaxies is consistent with values slightly larger than rjs = 0.5, which would correspond to a slight increase of the M/L-ratio with 
luminosity. The discontinuity at L = 10^'^ Lq is a numerical artefact, since we allow for two different scaling relations in the two 
luminosity intervals. We do not require that the two relations agree at the point of connection (the 1 a error intervals are shown as dashed 
lines). We included for comparison the single-power scaling for the combined sample in both panels in black {s oc L^'^'^, M/L ex. L^'^^, 
see Fig.[38j. 



best-fit slope for the c-L-scaling for given luminosity the 
mass of red galaxies is higher than for blue ones. Therefore 
the fiducial luminosity for the lens subsamples splits into 
two different fiducial masses. For massive galaxies, the 
concentration of red galaxies exceeds the concentration 
of their blue counterparts with same mass, but due to 
the steeper concentration-mass relation for blue galaxies, 
considering low-mass galaxies, this relation turns into its 
opposite. 

As we have seen in Section 14.21 in Figs. |3D] and [32] 
the scaling behavior of r2oo deviated from a power law 
when including galaxies with luminosities L < 10^" Lq. 
In principle there are two possibilities to explain this 
observation. Firstly one could assume a modification of the 
concentration-luminosity relation when moving to fainter 
galaxies, leading to increased values for r^oo, when not 
accounted for. Secondly the scaling behavior of the mass- 
to-light ratio could change, leading to higher masses and 
thus r2oo-values than expected with a single-power law. We 
will now further investigate both possibilities in more detail. 

To investigate a change in the concentration-luminosity 
relation, explicitly assuming that the M200 — L relation is 
described by a single-power law {M200 oc L for pure red 
and blue galaxies), we repeat our analysis and separately 
analyse the scaling behavior of the concentration parameter 



for galaxies brighter and fainter than 10^" Lq. We fix 
the fiducial virial radius r^oo including its scaling relation 
with luminosity (r2oo oc L^^^ for red and blue galaxies 
separately) and the fiducial concentration parameter c* 
and run maximum likelihood analyses, fitting independent 
slopes for the c-L relation for the brighter and for the 
fainter galaxies. Indeed we observe a different scaling 
behavior in both luminosity ranges. While for galaxies 
with L > 10^° Lq the scaling corresponds to what 
we obtained for the all luminosities sample, for the fainter 
galaxies a value of ?7c < — 1 is indicated. In order to exclude 
that we observe an effect of remaining unproperly treated 
multiple deflections, we repeat the analysis for galaxies in 
low density environments, obtaining the same result. This 
means, if we assume that the M2oo/i-ratio is described 
by a single-power law, the c-L relation comes significantly 
steeper for galaxies fainter than L = 10^° L®, leading 
to much higher values for the concentration. 

We further consider the second possibility, assuming a 
change of slope for the M2oo/I/-ratio (see earlier in this 
section) and thus for r2oo, as obs erved in Sec t ion I4.2T2] 
We assume the concentration of iDuffv et al.1 (|2008l ) to 
hold and fit the r2oo and its scaling index 7;r2oo fo'" bright 
and faint galaxies independently) and obtain a signifi- 
cantly shallower r2oo — L relation for galaxies fainter than 
Lj.1 = 10^" Lq, resulting in r2oo oc L" for red and 
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Figure 41. Scaling relation for the concentration parameter c with luminosity. The red contours show the scaling relation for red 
galaxies, the blue contours show the scaling relation for blue galaxies and the black contours for the combined galaxy sample. While 
for the combined and the red galaxy sample the result is consistent with a constant c-L-scaling, the best-fit value for the slope of the 
c-L-relation of blue galaxies is negative. However, the uncertainties are quite large and thus the confidence levels for r]'^"'^ and r;J?'"® 
overlap. The contours indicate the 68.3 per cent (solid), the 95.4 per cent (dashed) and the 97.3 per cent (dotted) confidence levels, the 
cross indicates the best-fit value. 



'"200'' oc LP'^^ for blue galaxies). This implies an increase 
in mass-to-light ratio for faint galaxies with decreasing 
luminosities which i s in consistence with the results of 
Ivan den Bosch et al.l (|2005h (see their fig. 3) . 

4.3.4 Companson to previous results 

iParker et all (|2007D analysed 22 deg^ of CFHTLS i'-band 
data, but without photometric redshift knowledge, choosing 
lens galaxies to have apparent magnitudes of 19 < i' < 22 
and source galaxies 22.5 < i' < 24. Assuming a trunca- 
tion rad ius of s = 185 ±30 kpc, based on the re- 
sults of iHoekstra. Yee fc GladdersI l|2004h . they found for 
a L^^ = 1.3 X 10^" h~^ Lq galaxy a velocity disper- 
sion of a* = 137 ± 11 km s~^ and a rest frame mass- 
to-light ratio of 173 ± 34 /i Mq/Lq. Taking the higher 
value for our fiducial luminosity into account, our values 
of a* = I31I2 km s~^ and our mass-to-light ratio of 
Mtotai.BBS /L — 1781^9 ^ -^0 / -^0 S'l'® slightly lower, however 



consistent. For a luminosity evolved mainly red lens s ample 
with luminosity Lr = 10^° /i"^ Lq Ivan Uitert et al.l (|201ll ) 
found a virial mass of Mvir = 7.2 ± 1.5 x lO" /i"^ Mq. If 
we convert our result for M200 from the NFW likelihood for 
a pure red sample to their reference luminosity we obtain 
a value of Mvir = 7.3 ± 0.5 x lO" Mp, , being in good 
agreement with the result of Ivan Uitert et all (|201ll ). 

4.3.5 Consistency of maximum likelihood and AS fit 
results 

We compare our results from the maximum likelihood anal- 
yses and the measurements of the excess surface density 
AE (see Tables [T] and (S)) and find in general a good agree- 
ment. For an L*-galaxy we obtain a velocity dispersion of 
a* = 135 ± 2 km s"^ from AE and a* = 13ll2 km s"^ 
from the maximum likelihood analysis for all galaxies. For 
blue galaxies we obtained values of o^^iue = 115 ± 3 km s~^ 
from AE and (jJiuo ~ II8I5 km s~^ from the likelihood 
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Figure 42. Excess surface surface mass density AS extracted from lens galaxies with —23 < M^i —24 calculated for simulated lens 
signals. The black dots show simulations including lenses with all magnitude, red dots only include lenses with M^/ —21, green dots 
only include lenses with M^i < —23, cyan dots additionally restrict the lens sample to maximum neighbour number of 4 galaxies within 
a projected distance of 720 kpc and finally the magenta plots only show those lenses without any neighbour within this projected 
distance. We see that in the case of a BBS profile (left panel) multiple deflection affect the signal amplitude only on larger scales, while 
assuming an NFW profile (right panel) the signal amplitude already at very low scale is biased high about 20 per cent. 



analysis. The consistency for the red lens sample with 
(7*3d = 162 ± 2 km s"^ from AE and a*^^ = 149t^ km s"^ 
from the maximum likelihood is only marginal. If we in- 
clude luminosity evolution, the results from the maximum 
likelihood analyses and from the AE-profiles in general 
also agree well. For all galaxies we obtained a velocity dis- 
persion of (T* — 150 ± 2 km s~^ from the AE-analysis and 
a* = 144;};2 km from the likelihood analysis, for blue 
galaxies we obtained (t^iuc = 123 ± 3 km s~^ from the AE- 
analysis and crj^iuo — 12(i^\ km s~^ from the maximum 
likelihood analysis. For red galaxies the values are again 
only marginally consistent with a*^^ — 173 ± 2 km s~^ from 
the AE-analysis and a*^^ — 161^3 km s~^ from the 
likelihood analysis. For the virial radius r2oo, the situa- 
tion is different, as in this case the results are in gen- 
eral only marginally consistent. The maximum likelihood 
analysis yields a value of rJoo = 133^2 kpc for the 
combined lens sample while the amplitude of the AE 
lens signal leads to a value of r2oo = 146 ± 2 kpc. 
Also for red (r^^^'^'* = 177 ± 3 ft-^ kpc from AE and 
rJoo"^ = I6OI4 kpc from maximum likelihood) and 
blue galaxies (''200'"'^ ~ ^20 ± 2 h^^ kpc from AE and 
rjoo'"" = 115^5 kpc from maximum likelihood) sepa- 
rately, the virial radii derived from the excess surface mass 
density are up to 10 per cent higher than in the maximum 
likelihood analysis. We obtain the same result when extract- 
ing the r2oo from a simulation-based AE-profile. 
In order to understand the origin of this apparent discrep- 
ancy, we perform 3D-LOS-projected lensing signal simula- 
tions (BBS and NFW), initially including all lenses. Incre- 
mentally we restrict the lens sample to M^' < —21, then we 



only allow lenses with Mr < —23. We then further restrict 
the lens galaxies to less than five neighbours within a pro- 
jected distance of 720 kpc and finally only consider iso- 
lated lenses without neighbours within the same projected 
distance. In this way we incrementally reduce the influence 
of multiple deflections (mainly by galaxies in the lens envi- 
ronment but also by other lens galax;ies in the line-of-sight). 
We use the simulated shape catalogues to calculate the cor- 
responding AE-profiles. For the analysis we focus on the lens 
luminosity bin with -24 < M^/ - 23 (see Fig. [IS]). The pro- 
files show for both models that on large scales the amplitude 
significantly decreases with decreasing importance of multi- 
ple deflections. This confirms the assumption that in this 
range the shear signal is significantly influenced by nearby 
galaxy haloes. Focussing on the BBS-based AE-profiles, the 
results for the various lens samples hardly differ in the in- 
ner regions and only deviate on larger scales. This implies 
that velocity dispersions derived from a single-halo-SIS-fit 
to the AE-profile on short scales instead of a full maximum 
likelihood analysis should be hardly biased. Surprisingly, the 
situation for the NFW-based profiles looks different. In this 
case multiple deflections lead to an increased signal ampli- 
tude also on short scales. I.e. virial radii derived from the 
AE in general are biased high about 20 per cent even on 
short scales, when not explicitly considering a lens sample 
in a low density environment. This can be explained, if we 
take into account, that the profile slopes for BBS and NFW 
behave differently. While the BBS profiles are isothermal on 
shorter scales and thus their slopes are independent of lu- 
minosity, this is not the case for NFW profiles. The brighter 
a lens, the larger its r2oo (and for similar concentration its 
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scale radius rs). For brighter galaxies, on short scales we 
basically only analyse regions within or a few scale radii, 
while for fainter galaxies the contribution of regions outside 
the Ts (where the profile slope becomes much steeper) signif- 
icantly increases. Therefore the slopes of the AE-amplitudes 
might increase with increasing contribution of fainter lenses. 
I.e. NFW halo parameters should not be purely or only care- 
fully extracted from AE profiles. 



5 SUMMARY AND CONCLUSIONS 

We performed a GGL analysis on 89 pointings of the 
CFHTLS-WIDE. We used the u*g'r'i'z' photometry to es- 
timate photometric redshifts with an accuracy of smaller 
than 0.04 (H-^) for objects with i' < 24.0. The background 
galaxy shap es were measu red and corrected with th e KSB-f- 
method (see lKaiseJll995l and lHoekstra et al.lll998l '). 



5.1 Weak lensing results from analyses on shorter 
scales 

We measured the tangential gravitational shear signal on 
scales from 25 kpc out to 700 kpc, examining 
all lens galaxies and subsequently red and blue galaxies 
separately. On average the measured velocity dispersions 
for red galaxies are about 25 per cent higher than for the 
average lens sample, while the velocity dispersions for blue 
galaxies are about 15 to 20 per cent lower. We also analysed 
different galaxy density environments. The tangential shear 
amplitude, especially on larger projected distances, shows 
higher values in denser environments, due to multiple 
deflection by nearby galaxy and by host haloes (see also 
lBrainerdl[2O10l ;i. 

We measured the excess surface mass density AE(/?) 
and compared our results for different luminosity in- 
terva ls with previous work, f inding partial consistence 
with iMandelbaum et al.l l|2006d ') . The results for the blue 
galaxy sample are in good agreement. The measurements 
for red galaxies show some differences in amplitude 
depending on t he luminosity bin. Our resu l ts ho wever 
agree weU with IMandelbaum. Seliak fc Hiratal (|2008l ). We 
further compared our measurements for red galaxies to 
Ivan Uitert et al.1 (|201ll ). finding good agreement for most 
investigated luminosity bins. 

We investigated the luminosity scaling relations for 
the velocity dispersion a, analysing AE in different lumi- 
nosity bins for our red, blue and combined lens samples 
on scales out to 100 h'^ kpc. We find a,^i oc L°-'^^^°-°^ 
for red and abiuc oc L ' ' for blue galaxies. For the 
same luminosities the red galaxies have a higher value 
of a than the blue ones. The combined galaxy sample is 
dominated by blue galaxies at low luminosities and red 
galaxies at high luminosities, which leads to a steeper 
relation a oci0-2»±0 02. 

Our results from fitting isolated haloes to AE(_R) are 
confirmed in a maximum likelihood analysis, thus prop- 
erly treating multiple deflections by independent lens 
galaxies. For an SIS profile out to maximum separations 



of 200 h'^ kpc, we obtain values for the a — L scaling 



of r]a 



0.31 



+0.02 



for all galaxies or 



0.27: 



0.27loQ4, when considering red and blue 



and -nT" 

SED galaxies separately. Therefore we have fixed the 
scaling relations to o" oc L" '' for the combined sample and 
(Trod/blue L°'^® for red and blue galaxy samples for the 
rest of our analysis. 

We analysed the scaling of the NFW-based halo parameters 
r2oo, M200 and V200, assum ing the concentration- mass- 
relation of iDuffv et al.1 (|2008l ). We find that for galaxies 
with rest frame luminosities L^i > W^'' Lq the 

scaling relations follow a power law. The power laws are 
raoo oc L^-^^^om ^^^^ ^ ^o.36±o.07 ^^j^^ 

galaxies. For the combined lens sample we find, analogously 
to the velocity dispersion a (and for the same reason), a 
steeper relation of r2oo oc 7;,0 39±0-04 rpj^^ virial mass of 
the combined lens sample scales as M200 oc /^i 2i±o.io 
However, for galaxies fainter than L = 10^" Lq the 
values of the NFW halo parameters remain significantly 
higher than suggested by a single-power law scaling 
relation. This either indicates a significantly slower de- 
crease of galaxy radii and masses or a significant increase 
of the concentration parameter c with decreasing luminosity. 

We compared the results for the velocity dispersions 
of o ur red len s samp le to kinematical me asurements from 
Gall azzi et al.l (j2006l ) and LRGs based on lEisenstein et al.l 
(;200ll ). We saw that 0.96 (Jetar ^ aZio ^ 1-07 cr^tar holds. 
This means that the halo velocity is between the maximal 
circular velocity found around ~ 0.5 Re and 90 per cent 
of this value which equals the circular velocity of galaxies 
at a few effective radii. Taken together this implies that 
haloes are isothermal from a few effective radii out to about 
100 kpc. Also, the values for the circular velocity V200 
of our blue galax y sa mple agree we l l with the results of 
iReves et al.l (j201ll ) and lDutton et al.1 (|2010l ). 

We find a break in the luminosity scaling of the M/L-ratio 
(red galaxies have a minimum M/L-ratio for a stellar mass 



of M,, 



3 - 4 X 10^" Mp,). This is in agreement 



with the result of iDutton et al.l (|2010l ). who combined weak 
lensing measurements with satellite kinematics. 



5.2 Weak lensing results from analyses on larger 
scales 

We performed maximum likelihood analyses as described 
by ISchneider fc Rixl l|l997h . We investigated two different 
halo profiles, an NFW profile (out to 400 kpc) and a 
truncated isothermal sphere (BBS, out to 2 Mpc). We 
assumed the scaling relations of to iGuzik fc Seljak (200 j) 

to scale the masses with luminosities and iDuffv etld] 

(|2008l ) for the mass-concentration relation, while further 
applying the scaling relations we previously extracted 
from the AE-profiles. For the maximum lens sample 
we find a virial radius of tIqq = 133^2 kpc and a 

concentration parameter of c* = 6.41q y with a viral mass 
of Afa'oo = 7.6t°o l x 10^^ h'^MQ. The virial radii are 
about 20 per cent higher for red galaxies and 15 per cent 
lower for blue galaxies. This results in virial masses about 
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60 per cent higher for the red galaxies and 35 per cent 
lower for blue galaxies. For the BBS profile we find a 
velocity dispersion of a* = 13f ita km and a truncation 
radius of s* = 184^^4 kpc. This corresponds to a 

total mass of Aft*„t,i,BBS = 2.32lg|8 x 10^" h"^ Mq. 
Both measurements are in agre ement with literature (see 
iHoekstra. Yee fc Gladdersl liooll). 

Considering red and blue galaxies separately, we find 
that red galaxies show significantly higher values for the 
velocity dispersion and the truncation radius than the blue 
galaxies. Blue galaxies have (for an L'-galaxy) basically 
the same concentration parameter as red galaxies but a 
significantly lower value for the virial radius. This leads 
to a significantly lower virial mass than for red galax- 
ies (about a factor of 2.5). The values for the best fitting 
halo parameters for an I/*-galaxy are summarised in Table[2l 

We analysed the properties of our lens galaxies as a 
function of their environment density. We split the galaxies 
into two halves relative to the median galaxy environment 
density. This results in a modestly underdense field galaxy 
dominated galaxy sample and a modestly overdense galaxy 
sample, where more galaxies reside in groups. We find that 
the velocity dispersion a of BBS haloes hardly depends on 
the environment, i.e. that the central galaxy halo profile 
mostly depends on the galaxy luminosity. The BBS trun- 
cation and NFW virial radii for red galaxies are however 
significantly larger for galaxies in denser environments. The 
increase of virial radii translates into 50 per cent higher 
virial masses for red galaxies in the moderately overdense 
relative to the moderately underdense sample. We do 
not discriminate between central and satellite galaxies in 
this work, i.e. our analysis is insensitive to the level at 
which satellite haloes are stripped and central galaxies 
grow in mass. We instead measure the average mass for a 
given luminosity in different environments. The increase of 
average halo masses is expected if one takes into account 
that the total mass-to-light ratio increases from galaxies 
to groups (and ga laxy clusters) (|van Uitcrt et al.. .20f I. . 
ISheldon et alll2009l ). 

We performed a BBS maximum likelihood analysis, fixing 
(T* and its luminosity scaling 770- to fit the truncation radius 
s* and its luminosity scaling relation a oc L^" . For the com- 
bined lens sample we find a value of t)s = 0.52^0%^ 

leading 

to a M/L-scaling of M/L oc L^-^^io l?. When we analysed 
the blue and red sample separately we found that the blue 
and the bright red sample have a M/L-scaling consistent 
with a constant M/L-ratio. For the faint (L < 10^" Lq) 
red galaxy sample we found M/L oc L^^ , with a minimum 
at ~ f0^° Lq (corresponding to M^tar ~ 10^° h~'^ Mq). 
The existence and location of this minim um is consistent 
wit h results from abun dance matching (see lGuo et al.llioiol 
and iDutton et aLll201(]| ') . 

We analysed the concentration-luminosity relation as 
a function of the combined lens sample and for red and blue 
galaxies separately. We find a scaling of c oc L ' -o-n 
for th e lens sample as a whole, consistent with iDuffv et al.l 
l|2008f ). The scaling of red galaxies is similar, but slightly 



flatter (c oc L 



), showing hardly any evolu- 



tion with luminosity or mass. For blue galaxies we flnd 

-0 34+°-^* 

c OC L ■ -0.26 , The best-fit value would imply a significant 
increase of the concentration with decreasing luminosity. 
This would mean, that at the high luminosity end of our 
lens sample the concentration of red galaxies exceeds the 
concentration of blue galaxies of same luminosity while for 
their faint counterparts the opposite is the case. However, 
the uncertainties for the measured slopes for blue galaxies 
are quite large. Therefore the values for ■ql^'^ and ?7c'"° still 
are consistent with each other within the uncertainties. 



For galaxies fainter than L — 10^'^ Lq we mea- 

sured higher values for r2oo from AE than expected for 
a single-power scaling with luminosity. We investigated 
two possibilities to fix this. On one hand we analysed the 
implications of a modified conc e ntrat ion- luminosity-relation 
in comparison to iDuffv et al.l l|2008l ). The analysis of the 
modified c-L scaling relation yields comparable values for 



bright 



as for the case where we considered all luminosities. 



For the faint luminosity range a slope of ry/'" < — 1 is sug- 
gested, leading to significantly higher concentration values. 
On the other hand we assumed a double-power-law scaling 
relation for the r2oo with luminosity. When allowing two 
different slopes for the r2oo-L-relation for L > fO^" Lq 
and L < h^^ Lq, the resulting boost of the fainter 

lenses' AE profiles makes an additional modification of the 
c-L scaling unnecessary and also gives better agreement to 
the observed AE-profiles for faint galaxies. The ch ange of 
slope in M200 /L is also predicted by ivan den Bosch et al.l 
(|2005l ). 



5.3 Comparison with 3D-LOS-projected lensing 
signal simulations 

Finally we created two sets of synthetic lensing observa- 
tions. The locations, the redshifts, the SED types and the 
luminosities of the lens and source galaxies were chosen to 
be identical to the true observations. The lensing properties 
of the lenses were described with either BBS or NFW 
profiles with parameters obtained from the maximum 
likelihood analyses regarding the Z/*-halo parameters and 
the luminosity scaling relations. 

We analysed the simulated observations in the same way as 
the true observations and measured -yt- and AE-profiles as 
a function of SED type and absolute luminosity and com- 
pared the results of both simulations with the observational 
data. The AE-profiles derived from the simulated obser- 
vations mostly agree well with the observed profiles. More 
precisely if the NFW simulation is based on a single-power 
law scalin g for r2oo and assu ming the concentration-mass 
relation of " 



Duffy et al.l (|2008f), it predicts a too low lensing 



signal for the fainter luminosity bins. This can be fixed by 
either assuming a steeper concentration-luminosity scaling 
or a shallower scaling relation for the r2oo of faint galaxies. 
The NFW simulation based on a double-power law scaling 
relation for the r2oo (^'200 

L^/^ for L > f0^° h-^ Lq, 
r'^lt oc L° and r^i^ « L"'^^ for L < 10'° ft'^ Lq) indeed 
leads to predictions agreeing significantly better with the 
observations. 
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APPENDIX A: SYSTEMATICS 

In this section we verify the integrity of our weak lensing 
analysis. We demonstrate that there is no weak lensing sig- 
nal if we replace foreground galaxies by stellar objects or 
random points and that our measured lensing signal ampli- 
tudes do not show a significant dependence on source mag- 
nitude, S/N-ratio or size. We also show, that misassignment 
of foreground objects to the background does not introduce 
a bias of more than ~ 1 per cent to the measurements of 
velocity dispersions. The contribution of faint not detected 
galaxies in the neighbourhood of more massive galaxies does 
not significantly affect the shear measurement of the massive 
galaxies and can be neglected. 



Al Gravitational shear estimates relative to stars 
and a random foreground distribution 

In absence of systematic errors the background galaxies 
must not show any systematic alignment relative to fore- 
ground stars or a random distribution of points. We there- 
fore take the stars used for our KSB PSF anisotropy cor- 
rection. These are 274 589 stars with magnitudes of about 
18 < i' < 22 and a SExtractor star classification of greater 
than 0.96. We treat this stellar sample as foreground sam- 
ple and analyse the lensing signal for our complete and un- 
affected background sample. As Fig. lAll (left panel) shows, 
both the tangential shear -yt and the cross-shear 7c are con- 
sistent with zero. In the second step we take all stars from 
our photometric catalogue with magnitudes 18 < i' < 22 
and a SExtractor star classification of greater than 0.96 
leading to total 'foreground' sample of 471 066 objects. Also 
in this case E-modes and B-modes are consistent with zero 
as can be seen in the right panel of Fig. lAll 
We then dice random points over the full CFHTLS area and 
then only consider random points outside the masked area. 
In this way the randomly positioned 'foreground' objects 
sample the same area as the original (photometric) fore- 
ground catalogue. We again keep the original background 
shape catalogue and measure E- and B-modes consistenct 
with zero for both quantities on all scales (see left panel 
of Fig. IA2[I . Finally we take the original photometric fore- 
ground sample and modify our background shape catalogue. 
We keep the original background positions but randomly re- 
assign the measured shears by permutation, conserving the 
ellipticity distribution. In this case the tangential shear and 
the cross-shear are in agreement with zero as well, see right 
panel of Fig. El 



A2 Estimate of the signal dilution from 

contamination of the background shape 
samples with foreground galaxies 

The lensing signal can be diluted if galaxies which are phys- 
ically related to the lenses are accidentally assigned as back- 
ground objects (e.g. due to photometric redshift mismatch). 
To quantify this effect we measure the density of back- 
ground objects around foreground objects within a radius of 
~ 2 h-^ Mpc as a function of t he projected separation (see 
also iMandelbaum et al.l l2005al ) . On short scales the num- 
ber density of background objects is slightly increased (see 
Fig. IA3|I by up to 3 per cent. Therefore we either have in- 
deed a dilution of the lensing signal by 1 to 3 per cent on 
scales 500 kpc or the increase in background galaxies 
is caused by the magnification effect. Using that 74 cx we 
conclude that in the worst case the velocity dispersions can 
be systematically underestimated by ~ 1.5 per cent on very 
short scales and less than ~ 1 per cent on larger scales due 
to contamination of background by foreground objects. 

A3 Signal contribution of undetected low mass 
nearby galaxies 

We now measure the contribution of low-mass objects to the 
lensing signal of the central galaxy. For larger redshifts these 
galaxies cannot be detected anymore and their halo mass is 
assigned to the central halo. To estimate the amount we 
run simulations (see Appendix [B] for details on the simu- 
lations). In the first setup we only consider bright galaxies 
(Mr' < —21) at lower redshifts (zphot ^ 0.5) with faint 
companions (M^/ > —21, maximal separation 720 Mpc 
in projected distance and 0.15 in redshift) but ignore the 
contribution of the companions in the lensing signal simu- 
lation. In the second step we also consider the lensing con- 
tribution of the companions. We then calculate the excess 
surface mass density AE for both cases. We see that the 
contribution of the low-luminosity companions is negligible 
(see Fig. lA4l red points for the bright galaxy signal vs. black 
points for the combined signal). For comparison we also in- 
vestigate the complementary lens sample, only comprising 
lenses with M^i > —21 (see Fig. IA4I blue points for the 
faint galaxy signal vs. black points for the combined signal). 
We conclude that for faint galaxies AE on large scales is 
dominated by mass associated with brighter galaxies in the 
neighbourhood. 

A4 Scaling of the lensing signal with redshift 

We now analyse the scaling behavior of the gravitational 
lens signal with redshift. We consider all source galaxies with 
0.6 Zphot ^ 2.0 and divide them into four different redshift 
bins of same size. We select lenses with of 0.05 < Zphot ^ 0.5 
to ensure the equality of the lens sample for independent 
background bins and avoid differences in the examined red- 
shift range due to our dynamical lens-source assignment. We 
then measure both the tangential shear signal and the ex- 
cess surface mass density and observe good agreement for 
all four subsamples (see upper left panel in Fig. IA5p . The 
velocity dispersions a for all four redshift bins, fitted within 
a projected separation of 200 kpc (see also Section [4. If) , 
are in agreement. The results are summarised in Table IXT] 
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Figure Al. Tangential gravitational shear (black symbols) and cross-shear (grey symbols) relative to a stellar foreground sample. Both 
the E-modes and B-modes are consistent with zero. The left panel shows the stellar selection performed with magnitude cut, SExtractor 
classification and the KSB-pipclinc, the right panel shows the stellar selection from the original photometric catalogues by selecting via 
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Figure A2. Tangential alignment of galaxies from the maximum background catalogue relative to random foreground positions (left 
panel) and tangential shear of ba<;kground with permuted shears relative to original foreground objects (right panel). The E-mode signal 
shown in black and the B-mode estimate in grey, are both consistent with zero. 



A5 Shear calibration biases for small and 
low-S/N objects 

We now use the complete foreground lens sample again. We 
split our background sources each into two subsamples of 
same number in terms of following properties: 
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Figure A4. Contribution of bright and faint galaxies to the excess surface mass density AS. The left panel shows the BBS simulations, 
the right panel shows the NFW simulations. The red triangles show the simulated signal contribution only for lenses with M^; < —21, 
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Figure A5. Systematic check. We split our samples into diflferent subsamples and compare the results. In the upper left panel we measure 
the excess surface mass density AS in for lenses with 0.05 < ^ 0.50 in four redshift bins. In the other three panels we consider the 
complete foreground sample. The upper right panel shows AS for bright (red symbols) and faint (blue symbols) sources, the lower left 
panel shows AS for large (red symbols) and small (blue symbols) sources and the lower right panel shows AS for sources with high (red 
symbols) and low (blue symbols) signal-to-noise. The measurements are in good agreement. 



(i) apparent magnitude, defining a bright and faint source 
sample, 

(ii) object size, defining a large and a small object sample, 
(ill) signal-to-noise ratio, defining a sample with high and 

low photometric signal-to-noise. 

We measure the gravitational shear and AE for all subsam- 
ples and obtain an estimate of the velocity dispersion u from 
an SIS fit to the shear. For the bright and for the faint 



background samples we measure a velocity dispersion of 
a — 115 ± 2 km and a = 119 ± 1 km s~^, respectively. 
For the large and for the small source samples we measure 
cr = 117 ± 2 km s~^ and a = 116 ± 2 km s~^, respectively, 
and for the high and the low signal-to-noise samples we 
measure velocity dispersions of cr = 115 ± 2 km s~^ and 
cr = 119 ± 1 km s~^. The AE estimates are shown in FiB:. IA5l 
to agree well with each other, see also Table IXT] 
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Table Al. 



name 






Background selection 


{Mr') 






iVfg 




N ■ 

^ " pair 


a [km s 


-^'bgl 


0.05 < ^fg 


s: 0.5 


0.60 ^ Zbg s; 0.73 


-21.0 


0.29 


0.67 


277 219 


292 900 


17 230 313 


121 ± 5 




0.05 < Zfg 


< 0.5 


0.73 ^ Zbg ^ 0.85 


-21.0 


0.32 


0.78 


277 577 


283 754 


19 174 161 


124 ±4 




0.05 < Zfg 


^ 0.5 


0.85 ^ Zbg ^ 1.01 


-21.0 


0.32 


0.92 


277 513 


273 113 


18 325 822 


118 ± 3 


Lbg4 


0.05 < Zfg 


^ 0.5 


1.01 ^ Zbg ^ 2.00 


-21.0 


0.32 


1.20 


277 598 


260 978 


17 531 415 


124 ±3 


-^bgbright 


0.05 < Zfg 


^ 1.0 


r' ^ 23.06 


-21.0 


0.38 


0.81 


1 711 502 


1 035 270 


67 033 496 


115 ± 2 


-^bgfaint 


0.05 < Zfg 


^ 1.0 


r' > 23.06 


-21.0 


0.48 


1.04 


1 811 810 


626 348 


67 025 398 


119 ± 1 


-^rhbig 


0.05 < Zfg 


^ 1.0 


rn =S 3.17 


-21.0 


0.40 


0.87 


1 772 769 


890 184 


67 031 724 


117 ±2 


-^rhsmall 


0.05 < Zfg 


^ 1.0 


Th > 3.17 


-21.0 


0.42 


0.92 


1 791 460 


771 434 


67 027 170 


116 ± 2 


-^snhigh 


0.05 < Zfg 


^ 1.0 


S/N 5s 9.34 


-21.0 


0.37 


0.80 


1 695 097 


1 037 322 


67 029 067 


115 ± 2 


-^snlow 


0.05 < Zfg 


^ 1.0 


S/N < 9.34 


-21.0 


0.48 


1.04 


1 819 346 


624 296 


67 029 827 


119 ± 1 



Results for the systematics check by splitting the background sample (redshift, apparent luminosity, object size and S/N-ratio) and 
analysing the subsamples with respect to the same foreground sample. The table shows effective lens luminosities and 
foreground/background-redshifts, number of lenses, sources and lens-source pairs and the velocity dispersions a from the SIS-fit to the 

inner region (R ^ 200 h^^ kpc.) 



APPENDIX B: 3D-LOS-PROJECTED LENSING 
SIGNAL SIMULATIONS 



The 3D-LOS-projected lensing signal simulations are based 
on our maximum likelihood results for the best-fit values of 
the halo parameters and their scaling relation with galaxy 
luminosity. 

For these shape-noise-free simulations we keep posi- 
tions, luminosities and redshifts from our previously defined 
lens and source samples and calculate the tangential shear 
signal 7t and the excess surface mass density AE induced 
by the lens sample to the source sample, in analogy to the 
maximum likelihood analysis, but with fixed values for the 
halo parameters. We use two different fiducial parameter 
sets for red and blue lens galaxies in our simulations. This 
takes the higher red galaxy mass for a given luminosity into 
account. Firstly assuming a truncated isothermal sphere 
(BBS) (see chapters l2.2.2l and l4.3.ip . we adopt our measured 
values for the velocity dispersion a*^^ =149 km s~^ and the 
truncation radius s*^^ = 337 kpc for red galaxies and 
ffbiue = 118 km s~^ and Sbiuc = 84 kpc for blue galax- 
ies. Secondly assuming an NFW profile (see Sections 12.2.31 
and I4.3.2|) . we adopt a concentration parameter c*ed ~ 6.4 
and a virial radius rJoo.rcd = 160 kpc for red galaxies 
and Cbiue ~ 7.0 and r-joo.biuc = US kpc for blue type 
galaxies. Since the simulations with a single-power law fit 
for the r2oo — L scaling relation lead to an underestimation 
of the excess surface mass density AE for the faintest 
luminosity bins based on the NFW profile, we assume a 
double-power law scaling relation for the r2oo with a change 
in slope at L = h'^ Lq (77^,00, bHght = <oo,bHght = 1/3, 
»?r2oo,faint = 0, VrlZ Mnt = 0.21). In addition we assume 
a change in slope for red galaxies for the s — L scaling 
behavior when considering a BBS profile (»7s^fiint ~ 0). 
We then calculate the gravitational shears for the galaxies 
expected from GGL, set them equal to the 'measured' 
ellipticities and rerun our analysis based on these artificial 
catalogues. 



Bl Predicted tangential alignment 

We repeat the analysis as described in Section [4. II and mea- 
sure the tangential shear 7t-profile for our simulated shear 
catalogue s. The qualitativ e results are in agreement to the 
results of lBrainerdI (|2010h . as on short scales we see an ex- 
cellent agreement of the profile to an SIS or BBS profile (the 
difference on scales significantly shorter than the truncation 
radius is negligible). On larger scales we observe a higher 
signal than expected from an SIS, due to multiply defiected 
sources by nearby galaxy haloes (see Figs. IBl1 and lB2] ). We 
observe that the reliability of the estimate for the velocity 
dispersion cr depends on the luminosity range of the inves- 
tigated lens sample. While for the brighter luminosity bins 
the results of the simulated tangential shear analyses agree 
quite well with the observations, for the fainter luminosity 
bins the simulated signal amplitude is slightly lower than 
the observed one. The reasons for this behavior can either 
be caused by not considered host haloes such as group or 
cluster haloes which are not sufficiently taken into account 
by only considering galaxy haloes. Another reason can be 
that the assumed scaling relations which are assumed to be 
power laws in luminosity are not universally valid but only 
a justified approximation depending on the investigated lu- 
minosity range. 



B2 Predicted versus measured excess surface 
mass density 

After the tangential shear profile analysis of the gravita- 
tional shear simulations we extend our investigation and cal- 
culate the excess surface mass density AS-profile (see Sec- 
tion based on our simulations. The results confirm the 
results from the tangential shear analysis: for bright galax- 
ies the signals from observational and simulated data agree 
very well down to a lens rest frame magnitude of M^/ ~ —20. 
For fainter (particularly BBS) lenses the amplitudes of the 
simulated sample starts to slightly underestimate the obser- 
vational signal increasing with decreasing luminosity, espe- 
cially on short scales. This effect appears to be weaker on 
scales larger than ~ 200 h~^ kpc, see also Fig. 1371 In general 
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Figure Bl. Tangential shear signal for simulated sample, BBS (left) and NFW (right) for galaxies with —24 ^5 M^i < —20. The colour 
code is the same as before, black stands for all galaxies, red for red galaxies and blue for blue galaxies. 
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Figure B2. Tangential shear signal for simulated samples, BBS (left) and NFW (right), looking at four individual one-magnitude 
luminosity bins with —24 ^ M^i < —20. 



the simulated BBS profiles fit tlie observational data better 
tfian the NFW profiles, when assuming single-power law fits 
for the scaling relations (not shown). However, when assum- 
ing double-power law fits for the scaling of the virial radius 
r2oo and the truncation radius s the NFW profile provides 
the better fit. 

If red and blue galaxies are considered separately, the simu- 
lations do not appear to fit equally well, in particular in the 



fainter luminosity bins. While the simulations of the blue 
and bright red (M^/ < —20) galaxies describe the observa- 
tional results fairly well, this is not the case for the fainter 
red galaxies, where the red lens simulations seem to par- 
tially underestimate the observed AE-profiles (see Figs. IB3I 
and|B4j. 

In the following we further analyse the simulated excess sur- 
face mass density as a function of the environment den- 
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Figure B3. Excess surface mass density AS for different luminosity bins only considering red galaxies. Black circles are the observational 
data points, red triangles come from the BBS simulation and green squares from the NFW simulation. 



sity. We split our lens sample as described in Section 13.91 
In the highest density environment (see Figs. IB5I and IB6|) . 
the simulations and the observational data match well for 
all luminosity bins. The profile slope for the lowest lumi- 
nosity intervals is rather flat, i.e. the dominating part of 
the lensing signal on scales larger than 200 kpc is not 
caused by galaxies from this specific luminosity bin but from 
brighter galaxies in the environment. The uncertainties for 
the faint low density environment lenses are quite large due 



to low number statistics. This makes it difficult to derive 
decisive conclusions. Nonetheless for the faintest luminosity 
bins the lensing signal is not inconsistent on scales shorter 
than R = 200 h'^ kpc (see Figs. [B7] and |B8ll. Taking into 
account the simple model we use to describe the mass asso- 
ciated to the galaxies, it is remarkable to see that the both 
simulations (based on the BBS and the NFW profile) agree 
very well with the observations over a wide scale of projected 
separations, galaxy luminosities and environment. Down to 
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Figure B4. Excess surface mass density AS for different luminosity bins only considering blue galaxies. Black circles are the observational 
data points, red symbols come from the BBS simulation and green squares from the NFW simulation. 



an absolute magnitude of M^i ~ —17 BBS and NFW sim- 
ulations mostly agree with each which makes it impossible 
to decide whether the observational galaxy haloes follow a 
BBS or an NFW profile. Only in the fainter bins the BBS 
and NFW signal start to differ, due to the differently mod- 
elled scaling relations for galaxies fainter than 10^° Lq. 
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Figure B5. Excess surface mass density for the lens sample in high density environments. Black circles are the observational data points, 
red triangles come from the BBS simulation and green squares from the NFW simulation. 
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Figure B6. Excess surface mass density for the lens sample in very high density environments. Black circles are the observational data 
points, red triangles come from the BBS simulation and green squares from the NFW simulation. 
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Figure B7. Excess surface mass density for tlie lens sample in low density environments. Black circles are the observational data points, 
red triangles come from the BBS simulation and green squares from the NFW simulation. 
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Figure B8. Excess surface mass density for the lens sample in very low density environments. Black circles are the observational data 
points, red triangles come from the BBS simulation and green squares from the NFW simulation. 



